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The Distance to SN 1999em in NGC 1637 from the Expanding Photosphere 

Method 

Douglas C. Leonard 1,2 , Alexei V. Filippenko 1 , Elinor L. Gates 3 , Weidong Li 1 , Ronald G. 
Eastman 4 ' 5 , Aaron J. Barth 6 , Schelte J. Bus 7 , Ryan Chornock 1 , Alison L. Coil 1 , Sabine Frink 8 , 
Carol A. Grady 9 , Alan W. Harris 10 , Matthew A. Malkan 11 , Thomas Matheson 1,12 , Andreas 

Quirrenbach 9 , and Richard R. Treffers 1 

ABSTRACT 

We present 30 optical spectra and 49 photometric epochs sampling the first 517 
days after discovery of supernova (SN) 1999em, and derive its distance through the 
expanding photosphere method (EPM). SN 1999em is shown to be a Type Il-plateau 
(II-P) event, with a photometric plateau lasting until about 100 days after explosion. 
We identify the dominant ions responsible for most of the absorption features seen in the 
optical portion of the spectrum during the plateau phase. Using the weakest unblended 
absorption features to estimate photospheric velocity, we find the distance to SN 1999em 
to be D = %.2±%% Mpc, with an explosion date of HJD 2,451,475.6 ± 1.4, or 5.3 ± 1.4 
days before discovery. We show that this distance estimate is about 10% closer than 
the distance that results if the strong Fe II AA4924, 5018, 5169 absorption features, 
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which have often been used in previous EPM studies, are used to estimate photospheric 
velocity. 

We examine potential sources of systematic error in EPM-derived distances, and 
find the most significant to result from uncertainty in the theoretical modeling of the 
flux distribution emitted by the SN photosphere (i.e., the "flux dilution factor"). We 
compare previously derived EPM distances to 8 SNe II in galaxies (or members of the 
same group) for which a recently revised Ccphcid distance exists from the HST Key 
Project and find -Dcepheids/-DEPM = 0.87 ± 0.06 (statistical); eliminating the 3 SNe II 
distances for which a Cepheid distance only exists to a group member, and not the 
host galaxy itself, yields -Dcepheids/-DEPM = 0.96 ±0.09. Additional direct comparisons, 
especially to spectroscopically and photometrically normal SNe II-P, will certainly help 
to produce a more robust comparison. 

Finally, we investigate the possible use of SNe II-P as standard candles and find 
that for 8 photometrically confirmed SNe II-P with previously derived EPM distances 
and SN 1999em, the mean plateau absolute brightness is My (plateau) = —16.4 ± 0.6 
mag, implying that distances good to ~ 30% (la) may be possible without the need 
for a complete EPM analysis. At M v (plateau) = -15.9 ± 0.2 mag, SN 1999em is 
somewhat fainter than the average SN II-P. The general consistency of absolute SNe 
II-P brightness during the plateau suggests that the standard candle assumption may 
allow SNe II-P to be viable cosmological beacons at z > 2. 

Subject headings: cosmology: observations — distance scale — galaxies: individual 
(NGC 1637) — supernovae: individual (SN 1999em) 

1. INTRODUCTION 

"Normal" core-collapse supernovae (SNe) are thought to result from isolated, massive stars 
(initial mass ^8 — 10 M & ) with thick hydrogen envelopes (generally several solar masses) intact 
at the time of the explosion. Their light curves show a distinct plateau (hence the moniker "SN 
II-P"), resulting from an enduring period (sometimes as long as 150 days) of nearly constant 
luminosity as the hydrogen recombination wave recedes through the envelope and slowly releases 
the energy deposited by the shock and by radioactive decay. It is to these events that the expanding 
photosphere method (EPM) of extragalactic distance determination is most accurately applied (e.g., 
Eastman, Schmidt, & Kirshner 1996, hereafter E96). 

SNe II have not traditionally been treated as standard candles. Through calibration with 
the EPM, however, they have provided direct measurements of distances independent of all the 
uncertain rungs of the extragalactic distance ladder. While EPM distances have been derived 
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to both SNe II-P and SNe II-L, 13 it is thought to be most securely applied to SNe II-P since 
the theoretical modeling of their flux distribution is much simpler and better understood. Thus 
far, the EPM has been applied to 18 SNe II, including 10 SNe II-P, spanning distances from 
0.049 Mpc (SN 1987A) to 180 Mpc (SN 1992am; Schmidt et al. 1994b). Using SNe II alone, a 
Cepheid-independent value of H Q = 73 ±7 km s" 1 Mpc" 1 has been derived (Schmidt et al. 1994a). 
Despite EPM's apparent success, however, there is considerable debate concerning possible sources 
of systematic error. In this paper, we derive an EPM distance to the most thoroughly observed SN 
II-P to date, SN 1999em, and critically evaluate possible sources of systematic error in the EPM 
technique. 

EPM is fundamentally a geometric technique, a variant of the Baade (1926) method used to 
determine the distance to variable stars: the linear radius of the expanding photosphere, R, is 
compared with the photosphere's angular size, 9, to derive the distance to the SN, D. Naturally, 
since all extragalactic SNe are unresolved during the "photospheric" phase (the plateau in a SN II- 
P), R and 9 must be derived rather than measured directly. 

In the original formulation of EPM (Kirshner & Kwan 1974), the theoretical angular size of 
the SN photosphere is approximated as 

» = §. a) 

where R the photosphere's radius and D the distance to the SN. If spherical geometry is assumed 
and the ejecta are in free expansion (i.e., deceleration due to gravity or swept-up interstellar material 
is negligible, and the radius of the progenitor is small compared with the size of the photosphere) 
equation (1) may be written as 

a _ V p \ 10 t(t - t ) 



D 



(2) 



where t> p hot is the velocity of material instantaneously at the photosphere at time t, and the explo- 
sion occurred at time t . Rewriting this equation in the suggestive form 

t = D (—) +t , (3) 



we see that plotted measurements of at least two epochs of #/v p hot against t should result in a line 
with slope D and y- intercept t Q . A particularly powerful feature of this method is that if the date 
of explosion (t Q ) is known (either from pre-discovery images or derived from the EPM technique 
itself if more than one epoch of data is available), then each measurement of 0/v p hot provides an 
independent estimate for D. To determine the distance to a SN, then, requires measurements of 
the theoretical angular size, 9, and the photospheric velocity, i> p hot- 



13 "L" for their linearly declining light curves, lacking a plateau; these events are generally believed to result from 
progenitors that have lost a substantial fraction of their hydrogen envelope prior to exploding (e.g., Filippenko 1997). 
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Assuming a spherically symmetric SN photosphere that radiates isotropically as a blackbody, 
conservation of flux demands (temporarily neglecting extinction) 



where B V {T C ) is the Planck function at color temperature T c , and f v is the flux received at Earth. 
Solving for 9, this becomes 



The observables thus required to derive an EPM distance are f p hot> fv, and T c . Traditionally, 
^phot is estimated by measuring the blueshift of the P-Cygni absorption troughs of weak lines seen 
in the optical spectrum, typically Fe II AA4924, 5018, 5169 and occasionally Sc II A5527 and Sc II 
A5658 as well. The implicit assumption is that these lines are optically thin above the photosphere 
and that the density of the absorbing material is a steeply declining function of radius, with the 
densest material (which causes the greatest amount of absorption) being instantaneously at the 
location of the photosphere. Since strong, optically thick lines may form absorption minima at 
higher velocities (i.e., well above the photosphere), velocities derived from the weakest features 
are the best indicators of photospheric velocity. In principle, f v and T c could be determined with 
accurate spectrophotometry. In practice, however, spectrophotometry is generally not available, 
so broadband photometry (usually some subset of BVIJHK) is used. From several measurement 
epochs of Vphot, T c , and f v during the photospheric phase, then, the distance and time of explosion 
of the SN are derived. 

The current debate over possible sources of systematic error in the EPM technique generally 
focuses on three main assumptions: 

1. The photosphere radiates as a blackbody. On applying the Baade (1926) technique to measure 
the distance to variable stars, Walter Baade himself wrote: "...the only assumption we had 
to make [in formulating the technique] was that Cepheids radiate as a blackbody. Further 
experience has to show how far we can hold this assumption." 14 It was not until decades later 
that Wagoner (1982) first demonstrated that a Type II SN photosphere in fact radiates as 
a "dilute" blackbody due to the dominance of electron scattering over absorption processes. 
That is, the continuum spectrum that is ultimately released from the electron-scattering pho- 
tosphere (defined as the surface of last scattering, r e = 2/3) is produced in the deeper layer 
at which the radiation field thermalizes to the local gas temperature, known as the "thermal- 
ization depth" (see, e.g., E96). This smaller radius is the location of the last true absorption 
and reemission of the photons that make up the thermal continuum. The blackbody spec- 
trum ultimately released by the electron-scattering photosphere thus possesses a luminosity 



4tt R\B U (T C ) = 47T.D 2 /. 



(4) 




(5) 



Passage translated from the original German by M. Modjaz. 
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appropriate for a smaller radiating surface 
by the relation 

c 

where -Rthcrm is the radius at the thermalization depth, i? p hot is the photosphere's radius 
(i.e., the surface of last scattering), and ( is called the "distance correction factor," since its 
inclusion in equation (4) "corrects" derived distances such that 

-^actual C -^measured- (7) 



The amount the flux is "diluted is parameterized 

-Rtherm /n \ 

= -B > ( 6 ) 

-Kphot 



That is, distances derived without taking into account flux dilution will overestimate the 
distance by a factor £. In principle, ( could depend on many things, including the chemical 
composition (i.e., the metallicity) and density structure of the progenitor star, and the ex- 
pansion rate and luminosity of the SN explosion. Studies of theoretical models of realistic 
SN atmospheres with a wide range of properties, however, have demonstrated that £ is in 
fact a nearly one-dimensional function of color temperature, T c , with only a small density 
dependence at short wavelengths (E96; Montes & Wagoner 1995). 

From a study of 63 model atmospheres, E96 provide convenient analytic approximations 
for C as a function of color temperature for a SN II-P, determined for 4 broadband filter 
combinations, BV, BVI, VI, and JHK (the R band is generally not used, since it is domi- 
nated by the Ha line, which can vary significantly among SNe II-P). In the models of E96, 
( is shown to vary smoothly as a function of T c , with a la scatter among the models of 
0.11%, 0.05%, 0.08%, and 0.07% for the BV,BVI,VI, and JHK filter combinations, respec- 
tively. Recently, by taking into account the effects that telluric features have on observed 
photometric data, as well as a subtle correction to the BVRI filter functions of Bessell [1990] 
to make them appropriate for use with photon, and not energy, distributions, Hamuy et al. 
(2001) has provided improved estimates of C(^c) ; which we adopt in this paper. 15 With BVI 
photometry of SN 1999em, we shall be able to test the self-consistency of the flux-dilution fac- 
tors by comparing the distances derived using the three filter combinations BV, BVI, and VI 
individually. 

The dilution factors of E96 have been used to determine distances to 18 SNe II (Schmidt, 
Kirshner, & Eastman 1992; Schmidt et al. 1994a, b; see also Clocchiatti et al. 1995). A 
first-order check on the precision of the predicted shape of the function relating ( and T c 
is whether multiple observations at different epochs, and therefore different T c and £, are 
consistent. As discussed by Leonard et al. (2001, hereafter L01), previous EPM studies of 



15 Hamuy et al. (2001) also derive an EPM distance to SN 1999em. We note that other than adopting the improved 
dilution factors our analysis here is completely independent of this work, which was posted on astro-ph after this paper 
was submitted; we do, however, make a brief comparison of the two derived distances, and comment on differences 
in the techniques used, in § 4.1. 
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SNe that had a broad temporal range of observational epochs have yielded consistent distances 
(i.e., a relatively straight line when 6/v is plotted against t [Equation (3)]). Empirical evidence 
therefore suggests that the shape of the theoretically derived function relating £ and T c must 
be nearly correct. However, differences of opinion remain concerning the overall level of the 
dilution factor, with discrepancies in £ as high as 60% claimed by some groups (e.g., Baron et 
al. 1995). In all, the uncertainty in the dilution factor remains the main point of contention 
when applying EPM to derive cosmological parameters. 

2. The photosphere is spherically symmetric. Recent polarization studies of young core-collapse 
SNe (see Wheeler [2000] for a comprehensive list of polarimetric observations of SNe; in 
general, higher polarization implies greater departure from a sphere), the possible association 
of SNe with some gamma-ray bursts (e.g., Bloom et al. 1999), and the high velocity of 
pulsars (Cordes & Chernoff 1998), all implicate a fundamentally asymmetric explosion for 
core-collapse SNe. From the polarization studies an interesting trend has also emerged: the 
degree of asphericity at early times appears to be inversely correlated with the amount of 
envelope material intact at the time of the explosion. That is, SNe that have lost a significant 
amount of their hydrogen envelope prior to exploding are found to be more highly polarized 
at early times than those that have hydrogen envelopes substantially intact (Wheeler 2000; 
see also Wang et al. 2001). The addition of envelope material evidently serves to dampen 
the asymmetry at early times, in agreement with the theoretical expectation that asymmetric 
explosions should turn spherical as they expand through the envelope (e.g., Chevalier & Soker 
1989; Hoflich, Wheeler, & Wang 1999). 

The most thorough investigation of the sphericity assumption for a SN II-P is that by L01 for 
SN 1999em itself. From multiple spectropolarimetric epochs spanning the first 163 days after 
discovery, L01 find a very low, but perhaps slowly increasing polarization for SN 1999em. 
Spectropolarimetric observations of other SNe II-P also show weak evidence for intrinsic 
polarization, suggesting that departures from spherical symmetry may not be substantial at 
early times for this class of objects (Leonard & Filippenko 2001). Although the translation 
from observed polarization to degree of asphericity is obfuscated by uncertainties in the 
interstellar polarization and SN viewing orientation, when the polarization of SN 1999em is 
interpreted in terms of the oblate, electron-scattering models of Hoflich (1991), a lower bound 
on the asphericity of 7% during the middle of the plateau results. If the asphericity implied by 
the spectropolarimetry results in directionally-dependent flux, it would clearly lead to distance 
errors in the EPM analysis for SN 1999em (see, e.g., Wagoner, 1991). We caution, however, 
that although global asphericity is the preferred interpretation of the observed polarization, 
there are other mechanisms that can produce SN polarization, including scattering by dust 
(e.g., Wang & Wheeler 1996) and asymmetrically distributed radioactive material within the 
SN envelope (e.g., Hoflich 1995; Hoflich, Khokhlov, & Wang 2001). Thus, while the question 
is far from settled, available polarimetric evidence is consistent with small asphericity during 
the period of interest for EPM. 
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Further support for early-time SN sphericity comes from the apparent success of the EPM 
technique itself. By quantifying the amount of systematic scatter in the EPM Hubble diagram, 
L01 find that the agreement between previous EPM distance measurements of SNe II and the 
distances to the host galaxies predicted by a linear Hubble law restrict mean SN II asphericity 
to values less than 30% (3<r) during the photospheric epoch. From the available evidence, 
then, there is growing support for the claim that asphericity does not significantly hamper 
the cosmological use of SNe II whose progenitors have massive hydrogen envelopes intact at 
the time of explosion. 

3. The photospheric velocity is accurately determined by measuring the blueshift of weak P-Cygni 
absorption lines. In all previous EPM studies, photospheric velocity has been estimated by 
measuring the blueshift of the flux minima in the P-Cygni absorption profiles of lines in the 
optical spectrum (e.g., Kirshner & Kwan 1974). The main assumptions here are that the lines 
in the SN envelope are optically thin above the photosphere, and that density is a steeply 
declining function of radius. We also note that due to the multiple scatters that occur to 
photons as they travel from the location of the thermalization depth to the surface of last 
scattering (i.e., r e = 2/3, the optical photosphere), line features are not expected to be formed 
significantly below the optical photosphere (i.e., Hatano et al. 1999; cf., Hamuy et al. [2001]). 
With these ideas in mind, then, it is clear that the greatest amount of absorption will occur 
just above the photosphere, and that velocities measured with weak lines should be good 
photospheric velocity indicators. 

The very early-time (less than 10 days after core collapse) spectra of SN 1987A were modeled 
in detail by Schmutz et al. (1990) and Eastman & Kirshner (1989). The hydrogen Balmer lines 
are shown by Eastman & Kirshner (1989) to yield velocities that become a progressively poorer 
representation of the photospheric velocity with time: they are within 10% of the photospheric 
velocity two days after collapse, but by ten days Ha, H/3, and H7 yield velocities 1.9, 1.5, and 
1.4 times the photospheric velocity, respectively. The explanation for this behavior comes 
from examining the density structure of the material overlying the photosphere. At early 
times, the atmosphere is highly ionized and the photosphere is situated in a region with a 
density that is a very steeply declining function of radius (e.g., Matzner & McKee 1999). 
As the atmosphere cools and starts to recombine the lines form in a much more extended 
region characterized by a flatter density gradient (Eastman & Kirshner 1989). The flatter 
the density gradient the more extended the line-forming region becomes, which produces the 
larger blueshift of the optically thick lines. At very early times, then, we expect all lines 
to converge on the photospheric velocity, but quite soon thereafter weaker features must 
be used. Both Schmutz et al. (1990) and Eastman & Kirshner (1989) conclude that the 
weak Fe II AA4924, 5018, 5169 lines provide reliable photospheric velocities during the early 
period studied for SN 1987A. Encouraged by these results, these features have been used 
as photospheric velocity indicators in most subsequent EPM distance determinations (e.g., 
Schmidt et al. 1992; Schmidt et al. 1994a, b; Clocchiatti et al. 1995), with the occasional 
use of Sc II A5527 and Sc II A5658 as well. Most of these EPM studies, however, used data 
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from significantly later evolutionary epochs than those used for SN 1987A. 

It is possible that the Fe II AA4924, 5018, 5169 features may themselves become optically 
thick above the photosphere as the temperature cools to the recombination temperature of 
hydrogen. Hatano et al. (1999) derive the expected local thermodynamic equilibrium (LTE) 
optical depths at the photosphere as a function of temperature for prominent absorption lines 
from each of the major ions thought to be responsible for features in a SN atmosphere with 
solar metallicity. The predicted depths nicely reflect the order of appearance of lines generally 
seen in SNe II-P as they cool: first hydrogen Balmer (which is shown to become optically 
thick below 12,000 K) and He I A5876 followed by Fe II, Si II, Sc II, and then other weaker 
metal lines. Using results derived in § 2 and § 3, Figure 1 illustrates the color-temperature (in 
this case, the blackbody temperatures determined by the BVI photometry, Tbvi) evolution 
of line strengths for SN 1999em. The line strengths and order of appearance are in good 
qualitative agreement with the predictions of Hatano et al. (1999); of particular interest, 
Hatano et al. (1999) predict Fe II A5018 to be optically thick above the photosphere at 
temperatures below T e g ~ 7500 K, suggesting that velocities derived using this line during 
the recombination phase may slightly overestimate photospheric velocity. By comparing the 
velocity derived from the Fe II AA4924, 5018, 5169 features with that inferred from weaker 
lines, we shall empirically test this for SN 1999em. 

With the preceding discussion in mind, we now explicitly include the effects of extinction (A v ) 
and flux dilution (Q in equation (4) and rewrite equation (5) as 



If accurate spectrophotometry were available, we could solve for 9 and T c by minimizing the quantity 
(f u — ( 2 2 ttB u [T c ]10~ oaa ' / ) over the observed frequency range. Since such data are seldom obtained 
we recast equation (8) in terms of broadband photometry. Following E96, we replace f v with the 
observed broadband magnitude (for brevity writing only BVI here) to arrive at 



where bsvi is the broadband magnitude of the Planck function integrated over the transmission 
function for the B, V,or I bandpass (plus an arbitrary constant). Armed with an estimate of the 
extinction and a set of photometric magnitudes, we solve for 6 and T c by minimizing the quantity 



U = cle\B v (T c )iv 



-0.4A„ 



(8) 




(9) 



e = 



J2i^BVi ~ A BVI + 51og# + 51og[C(T c )] - b BV i(T c )} 2 . 



(10) 



BVI 



With Vphot measured from spectral lines and 9 derived from photometry, the distance and explosion 
date of the SN are then determined through equation (3) by finding the best-fitting line to describe 
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14000 12000 10000 8000 6000 

Temperature (K) 

Fig. 1. — The strength of several optical absorption lines plotted against BVI color temperature 
for SN 1999em. Approximate spline fits to the observations are shown for clarity. Line strength 
is set equal to the measured line depth, defined as (f c — / m in)//c> where f m i n is the flux at the 
line's minimum and f c is the value of the flux in the surrounding regions interpolated over the 
absorption line (§ 2.2.3). The order of appearance and relative strengths of the lines are in good 
qualitative agreement with the theoretical predictions of Hatano et al. (1999). Note the presence 
of Si II A6355 absorption at high photospheric temperature (Tbvi > 10, 000K) and then again when 
the photosphere has cooled to below Tbvi ~ 6, 000K. 
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the data, either through the method of least squares or a more robust (i.e., less sensitive to outliers) 
estimator, such as the criterion of least absolute deviations (e.g., Press et al. 1992). 

SN 1999em was discovered on 1999 October 29.44 UT (UT dates are used throughout this 
paper) by Li (1999) at an unfiltered magnitude of m ~ 13.5 mag in the nearly face-on (i ~ 32°, from 
LEDA 16 ) SBc galaxy NGC 1637 (Fig. 2). Since an image of the same field on Oct. 20.45 showed 
nothing at the position of SN 1999em (limiting mag about 19.0; Li 1999), it was likely discovered 
shortly after explosion. It was quickly identified as a Type II event, with prominent P-Cygni 
features of hydrogen Balmer and He I A5876 and a blue continuum (Jha et al. 1999). Immediately 
after discovery a spectral, spectropolarimetric, and photometric campaign was initiated at Lick 
Observatory. Nearly daily optical photometry and spectroscopy of the young SN were obtained 
with the 0.8 m Katzman Automatic Imaging Telescope (KAIT; Treffers et al. 1997; Li et al. 
2000; A. V. Filippenko et al., in preparation) and the Lick 1-m Nickel reflector, respectively, 
providing unprecedented temporal coverage of this Type II event during the photospheric phase. 
The spectropolarimetric observations are discussed in detail by L01; here we focus on the spectral 
and photometric data. 

In this paper we derive an EPM distance to SN 1999em, and present 30 optical spectra and 
49 photometric epochs spanning the first 517 days of its development. We present the photometric 
data in § 2.1, and the total flux spectra in § 2.2, where we also identify and discuss the dominant ions 
responsible for the absorption features seen in the optical spectrum during the plateau phase. We 
consider the reddening of SN 1999em in § 2.3. We derive and compare the photospheric velocities 
inferred from the P Cygni absorption minima of several line features in § 2.4, and apply the EPM 
to SN 1999em in § 3. We discuss the results in § 4, and summarize our main conclusions in § 5. 

2. Reductions and Analysis 
2.1. Photometry 

All photometric data were collected using the KAIT, which is equipped with a 512 x 512 pixel 
Apogee CCD camera (AP7) located at the //8.17 Cassegrain focus, providing a field of view of 
6'.7 x 6'.7 with O'.'S per pixel. The "seeing", estimated by the full width at half maximum (FWHM) 
of stars on the CCD frame, generally ranged from 3" to 4'/5, and exposure times of 3 to 4 minutes 
in VRI and 5 to 7 minutes in UB were typical during the first 95 days after discovery while the SN 
was bright and on the plateau portion of its light curve. After SN 1999em dropped off the plateau, 
exposure times were increased to about 10 minutes. 

CCD frames were flatfielded using twilight-sky flats in the usual manner, and cosmic rays 
were removed interactively by visual inspection of each frame. A noticeable low-level pattern was 



http://www-obs.univ-lyonl.fr/leda/home_leda.html. 
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removed from each image by subtracting a dark frame scaled to the same exposure time; however, 
on some images, particularly the longer U and B exposures, a pattern still remained, most likely 
caused by temperature fluctuations of the thermoelectrically cooled CCD camera (Modjaz et al. 
2001). This was removed by interactively adding or subtracting an additional fraction (generally 
< 10%) of the dark frame, until the pattern disappeared. Photometric tests performed on frames 
before and after the additional dark subtraction indicate that the low-level pattern has a negligible 
effect on the photometry. Considerable fringing remained in /-band images that did not properly 
flatten, likely produced by the varying intensity of night-sky emission lines. This also has minimal 
impact on the photometry since both the SN and the comparison stars were much brighter than 
the background. 

Figure 2 shows a KAIT 5-band image of NGC 1637 taken on 1999 November 1. Of the 8 
"local standards" identified in the field of SN 1999em, the best were judged to be the three bright, 
isolated stars labeled A, B, and C, and they were used to measure the relative SN brightness 
on non-photometric nights. The absolute calibration of the field was accomplished on the two 
photometric nights of 2000 April 3 and September 25 by observing several fields of Landolt (1992) 
standards over a range of airmasses. The mean UBVRI magnitudes of the stars on the two nights 
are given in Table 1. The transformation coefficients to the standard Johnson-Cousins (Johnson et 
al. 1966 for UBV; Cousins 1981 for RI) systems were derived using the solutions from these nights 
and two additional photometric nights, 2000 September 5 and 27. We determined the instrumental 
magnitudes for the standards using aperture photometry with the IRAF 17 DAOPHOT package 
(Stetson 1987, 1991), which yielded color terms for the KAIT observations of the form 

U = u + 0.074(£/ - B) + Cu, 

B = 6 + 0.06(5 - V) + C B , 

V = v - 0.04(5 - V) + C v , (11) 

R = r - 0.08(V - R) + C R , 

I = i + 0M(V-I) +C7, 

where ubvri are the instrumental and UBVRI the standard Johnson-Cousins magnitudes. The 
terms Cu,Cb,Cv,Cr, and Cj are the differences between the zero-points of the instrumental 
and standard magnitudes, determined for each observation by measuring the offset between the 
instrumental and standard magnitudes and colors of the local standard stars. 

The location of SN 1999em at about 23" from the center of NGC 1637 (Figure 2) on the inner 
edge of a spiral arm gives it a rather complex background. In such situations, it is desirable to wait 
until the SN fades beyond detection to obtain a "template" observation of the host galaxy, which 
can then be subtracted from the earlier observations (Filippenko et al. 1986; Hamuy et al. 1994; 



17 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 
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Fig. 2. — i3-band image of NGC 1637 taken on 1999 November 1 with the Katzman Automatic 
Imaging Telescope. SN 1999em (SN) is measured to be 15'.'l west and if! 2 south of the galaxy 
nucleus (cf. Jha et al. 1999). The local standards listed in Table 1 are marked. 
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Table 1. Magnitudes of Local Standards 



Star 


U 


B 


V 


R 


I 


A 


14.97 


13.63 


12.42 


11.75 


11.17 


B 




16.75 


15.77 


15.18 


14.69 


C 


15.66 


15.74 


15.12 


14.76 


14.38 


D 


15.96 


16.21 


15.61 


15.25 


14.89 


E 


15.44 


14.16 


12.99 


12.37 


11.84 


F 


15.21 


14.99 


14.30 


13.94 


13.59 


G 




14.22 


13.50 


13.09 


12.72 


H 


13.85 


13.51 


12.79 


12.39 


12.04 
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Richmond et al. 1995), leaving a "clean" image of the SN with no galaxy contamination. However, 
nearly a year after discovery SN1999em still remained quite detectable in 10-min KAIT exposures, 
especially in the VRI bandpasses. Fortunately, during the period of main interest to this study 
(the first 100 days), SN 1999em was over 6 mag brighter than the galaxy background, reducing 
uncertainty due to improper background removal to at most the few percent level (e.g., Suntzeff et 
al. 1999). To obtain the most accurate photometric measurements without the benefits of galaxy 
subtraction, we used point-spread function (PSF) fitting by determining a PSF for each image with 
DAOPHOT using bright (but unsaturated) isolated stars in the field. We used only the inner core 
of SN 1999em and the local standards to fit the PSF in order to reduce errors that can be introduced 
when there is a strong gradient in the background (e.g., Schmidt et al. 1993). In practice, this core 
was generally set to be about the FWHM of a given image. 

One advantage of using the PSF fitting technique over simple aperture photometry is that it 
produces an image with the SN and local standard stars subtracted away, which allows a visual 
check on whether the brightness of the SN has been accurately measured. While the fitting radius 
of the SN and comparison stars was varied from night to night to match the seeing, the sky 
background of the SN and local standards was always set to an annulus with a radius of 20 — 26 
pixels (16'.'0 — 20 / /8) to maintain consistency throughout the observations. We then subtracted 
the mode of the sky background 18 to derive the instrumental magnitudes for the SN and local 
standards, and transformed them to the standard Johnson-Cousins system by taking the weighted 
mean of the values obtained using the three calibrator stars. Due to the faintness of the field stars 
in the U band, only star A was used to transform the ?7-band magnitudes. This limitation for U, 
coupled with the general difficulty of calibrating the U band using traditional CCDs with poor blue 
response (see, e.g., Suntzeff 2000), makes us regard the U photometry as systematically suspect at 
the 0.1 — 0.2 mag level. 

The results of our photometric observations are given in Table 2 and shown in Figures 3 and 
4. The reported uncertainties come from two sources. First, there is the photometric uncertainty 
reported by the error analysis package in DAOPHOT from the statistics of the SN and background 
region. Second, there is uncertainty in the transformation to the standard system. We estimate the 
transformation error by taking the standard deviation of the spread of the standard magnitudes 
obtained using each of the local standard stars. Since three values do not determine a robust stan- 
dard deviation, we opted to take a conservative estimate of the transformation error by measuring 
the standard deviation of the values obtained from all eight of the stars labeled in Figure 2. Due 
to the faintness of the local standards in the U band, the transformation error in that bandpass 
was generally obtained using fewer stars. The photometric and transformation errors were then 
added in quadrature to obtain the uncertainty reported in Table 2; in most cases, the total error 
was dominated by the uncertainty in the transformation. 



18 See Da Costa (1992) for a discussion of the advantages of using the mode of the background region rather than 
the mean or median. 
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Table 2. Photometric Observations of SN 1999em - Continued 



UT Date a Day b U (au) B (a B ) V (a v ) R (a R ) I (<tj) 
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a yyyy-mm-dd . 

b Day since discovery, 1999-10-29 UT (HJD 2,451,480.94). 
c Cloudy. 

d Poor seeing (^ 6"). 

e Poor guiding - images trailed. 

f Bright moon. 
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Fig. 3. — UBVRI light curves for SN 1999em from Table 2. For clarity, the magnitude scales for 
UBRI have been shifted by the amounts indicated. 
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Fig. 4. — U — B, B — V, V — R, and V — I color curves of SN 1999em. The magnitude scales for 
U — B, V — R, V — I have been shifted by the amounts indicated. 
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The VRI light curves of SN 1999em shown in Figure 3 indicate a clear plateau lasting until 
about 95 days from discovery. The plateau phase is characterized by a steeply declining L^-band 
brightness, a declining i?-band brightness, a slightly declining V-band brightness, a nearly constant 
-R-band brightness, and an /-band brightness that increases through the first 70 days and then 
slowly declines for the remainder of the plateau. The [/-band and B-band light curves decline by 
an average of 0.06 and 0.02 mag day -1 , respectively, during the plateau epoch after an initial quick 
rise to maximum on about day 2 after discovery. We measure the decline in the B band over the 
first 100 days after maximum to be fif^ = 2.1 mag, firmly establishing SN 1999em as a Type II-P 
event according to the definition of Patat et al. (1994), used to discriminate between SNe II-P 

(Pioo < 3 - 5 ma s) and SNe n " L (Pioo > 35 ma g)- 

The UBVRI photometric behavior of SN 1999em is broadly understood by considering the 
temporal evolution of the radius at the thermalization depth and color temperature shown in 
Figure 5 (using quantities subsequently derived by the EPM analysis in § 3). The surface area at 
the thermalization depth is seen to increase throughout most of the plateau as the overall expansion 
of the ejecta competes with the inwardly moving recombination front. If this were the only effect 
operating, we would expect to observe increasing brightness throughout the majority of the plateau 
phase. Countering this, however, is the decreasing temperature of the photosphere as the SN cools, 
which lowers the flux at all wavelengths. The effect of the decreasing photospheric temperature is 
least pronounced at longer wavelengths (i.e., the / band) since bandpasses significantly redward 
of the blackbody peak are more nearly in the Rayleigh-Jeans regime of the Planck function and 
thus only linearly dependent on temperature. As one moves to progressively bluer wavelengths 
the effect of the decreasing temperature becomes more pronounced, and the plateau brightness 
changes from generally increasing (/ band) to nearly constant (R band), to slightly declining (V 
band), to moderately declining (B band), to steeply declining (U band); significant line blanketing 
by metals further reduces the flux in the U and B bandpasses. The precipitous drop observed 
in all bandpasses after about day 95 signals the point at which the hydrogen recombination wave 
has receded completely through the massive hydrogen envelope and the transition to the nebular 
phase has begun. The full nebular phase, usually defined as the point at which the total electron- 
scattering optical depth to the center drops below unity, is not expected until after about day 250 
(Eastman et al. 1994). 

One feature of the F-band light curve (and possibly R as well) that warrants special comment 
is the distinct "double peak" evident in the early stages of development. The first peak occurs 
near day 4 after discovery (perhaps a bit later than the single peaks seen in U and B) and the 
second peak of nearly equal brightness occurs on day 13, with the dip between them reaching a 
local minimum on day 9 (0.07 mag fainter than the peak brightness). Similar photometric behavior 
can be discerned in the light curve presented by Schmidt et al. (1993) for the Type II-P SN 1990E, 
although in that case the second peak is brighter than the first, and occurs about two weeks after 
the first peak. In both SN 1999em and SN 1990E the rise to the second maximum coincides 
with the appearance of Fe II absorption lines in the spectrum (§ 2.2) and the point at which the 
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Fig. 5. — The temporal evolution of the radius of the radiating surface at the thermalization depth 
and BVI color temperature for SN 1999em during the photospheric phase (i.e., R = -Rthcrm = &(D 
with the theoretical angular size [#bvi]> dilution factor [Cbvi]> an d color temperature Pbvi] taken 
from Table 7, and a distance of D = Dbvi = 8.27 Mpc [Table 8] assumed). 
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photospheric color temperature has cooled down to near 10,000 K. This effectively marks the onset 
of the recombination phase. It is possible that the additional liberation of recombination energy 
slows the fall of the effective temperature enough to allow the effects of the increasing photospheric 
surface area to temporarily dominate the visual brightness of the SN. Additional examples of well- 
sampled early-time SNe II-P light curves are needed to see if this effect is generic. 

2.2. Spectra 

2.2.1. Observations, Reductions, and Identification of Line Features 

Table 3 lists the spectral observations of SN 1999em. With the exception of the single echelle 
spectrum, whose reduction details will be given elsewhere (D. C. Leonard et al., in preparation), all 
one-dimensional sky-subtracted spectra were extracted optimally (Home 1986) in the usual manner, 
generally with a width of ~ 10" along the slit, although this was modified based on seeing conditions 
(usually a width of ~ 4 times the seeing was employed) . Each spectrum was then wavelength and 
flux calibrated, as well as corrected for continuum atmospheric extinction and telluric absorption 
bands (Wade & Home 1988; Bessell 1999; Matheson et al. 2000). 

In general, the position angle of the slit was not aligned along the parallactic angle (the slit 
of the Lick 1-m spectrograph is fixed at 0°) so the spectral shape may suffer from differential light 
loss (Filippenko 1982). In addition, the shape of all spectra obtained with the Lick 1-m reflector 
is suspect, especially at wavelengths greater than 5800 A, since successive observations from the 
same night occasionally showed significant variations. 19 Multiple extractions of the data using 
apertures of varying width and different background regions failed to remedy the problem, and its 
cause remains unknown. Fortunately, this calibration uncertainty has minimal scientific impact 
here since it is the wavelengths of absorption features, not absolute fluxes, that are of primary 
interest. In any case, the effect was evident on only a few nights, and quantities derived using 
spectra obtained on these nights are in good agreement with those derived from spectra taken at 
other telescopes. 

The spectral evolution of SN 1999em during the first 517 days of its development is shown in 
Figures 6 and 7. The early spectra are characterized by a nearly featureless continuum with broad 
hydrogen Balmer and He I A5876 P-Cygni lines. As early as day 7, a hint of Fe II A5169 absorption 
is visible, and it becomes quite strong by day 11 along with Fe II A5018. During the next 90 days 
the strength of numerous metal lines continues to increase (see Figure 1). After the plateau ends, 
the spectrum becomes more emission dominated and marked by the presence of such forbidden 
lines as [Call] AA7291,7324 and [O I] AA6300, 6364. 



19 Normalizing the spectrum at the blue end, the flux level of the red end varied by as much as ±30% on the worst 
nights. 



Table 3. Journal of Spectroscopic Observations of SN 1999em 



Day a 


UT Date 


Tel. b 


Range c 


Res. d 


P.A. C 


Opt. P.A. f 


Air. s 


Flux Std. h 


See. 1 


Slit 


Exp. 


Observer (s) J 




(yyyy-mm-dd) 




(A) 


(A) 


(deg) 


(deg) 






(arcsec) 


(arcsec) 


(s) 




1.9 


1999-10-31 


LI 


3948-7449 


6 





150 


1.5 


HD19 


2.0 


2.9 


2x900 


EG 


2.9 


1999-11-01 


LI 


3948-7449 


7 





147 


1.5 


HD19 


2.2 


2.9 


3x900 


EG 


3.9 


1999-11-02 


LI 


3948-7449 


7 





154 


1.6 


HD19 


1.4 


2.9 


2x900 


EG 


4.0 


1999-11-02 


L3h 


3630-10100 


0.1 






1.3 




2.2 


2.0 


5400 


AQ, SF, CG 


4.9 


1999-11-03 


LI 


3948-7449 


6 





155 


1.5 


HD19 


2.4 


2.9 


3x900 


EG 


5.9 


1999-11-04 


LI 


3948-7449 


6 





159 


1.4 


HD19 


1.8 


2.9 


3x900 


EG 


7.0 


1999-11-05 


L3p 


3400-5500 


5 


170 


174 


1.4 


F34 


2.3 


2.0 


4x1200 


AF, DL, RG 








5000-7700 


6 


170 


174 


1.4 


HD84 


2.3 


2.0 


4x1200 


AF, DL, RC 








3400-5500 


5 


20 


24 


1.5 


F34 


2.3 


2.0 


4x1500 


AF, DL, RC 








5000-7700 


6 


20 


24 


1.5 


HD84 


2.3 


2.0 


4x1500 


AF, DL, RC 


7.9 


1999-11-06 


LI 


3948-7449 


6 





157 


1.4 


HD19 


2.1 


2.9 


3x900 


EG 


8.9 


1999-11-07 


LI 


3948-7449 


8 





158 


1.4 


HD19 


4.2 


2.9 


3x900 


EG 


10.2 


1999-11-08 


KII 


4000-10000 


14 


63 


61 


1.6 


BD17 


1 


0.7 


30+70+70 


AF, AC, AR 


11.2 


1999-11-09 


KII 


4000-10000 


14 


58 


61 


1.5 


BD17 


1 


0.7 


60+70 


AF, AC, AR 


20.9 


1999-11-19 


LI 


3948-7449 


7 





159 


1.4 


HD19 


3 


2.9 


3x900 


EG 


23.9 


1999-11-22 


LI 


3948-7449 


7 





162 


1.4 


HD19 


2.2 


2.9 


4x900 


EG 


25.9 


1999-11-24 


LI 


3948-7449 


7 





166 


1.4 


HD19 


2.7 


2.9 


4x900 


EG 


29.9 


1999-11-28 


LI 


3948-7449 


7 





164 


1.4 


HD19 


3.0 


2.9 


4x900 


EG 


33.9 


1999-12-02 


LI 


3948-7449 


7 





160 


1.4 


HD19 


3.5 


2.9 


4x900 


EG 


37.1 


1999-12-05 


KII 


4000-10000 


14 


61 


61 


1.0 


HD84 


2.2 


0.7 


330 


AF, PG 


39.9 


1999-12-08 


L3p 


4230-7010 


7 


16 


28 


1.5 


HD84 


3 


3.0 


4x1200 


DL, DA 








4230-7010 


7 


43 


45 


2.1 


HD84 


4 


3.0 


4x1200 


DL, DA 


43.9 


1999-12-12 


LI 


3948-7449 


7 





177 


1.4 


HD19 


2.6 


2.9 


3x900 


EG 








5550-8649 


7 





177 


1.4 


HD19 


2.6 


2.9 


3x900 


EG 


46.8 


1999-12-15 


LI 


3948-7449 


7 





177 


1.4 


HD19 


2.7 


2.9 


3x900 


EG 








5550-8649 


7 





177 


1.4 


HD19 


2.7 


2.9 


3x900 


EG 


48.8 


1999-12-17 


L3p 


4220-6882 


6 


161 


170 


1.4 


HD19 


2.3 


2.0 


4x1800 


DL, RC, AC 








4220-6882 


6 


24 


24 


1.5 


HD19 


2.5 


2.0 


4x1800 


DL, RC, AC 


75.9 


2000-01-13 


P 


3500-5500 


11 


90 


43 


1.7 


HD84 


2.1 


1.0 


3x120 


AH, SB 








5400-10000 


14 


90 


43 


1.7 


HD84 


2.1 


1.0 


3x120 


AH, SB 








3500-5500 


27 


90 


43 


1.7 


HD84 


2.1 


6.0 


2x60 


AH, SB 








5400-10000 


53 


90 


43 


1.7 


HD84 


2.1 


6.0 


2x60 


AH, SB 


94.8 


2000-02-01 


KII 


4510-7010 


5 


46 


141 


1.3 


G191 


0.9 


1.0 


120 


WB, DS, AS, WV, MR, 


95.8 


2000-02-02 


LI 


3948-7449 


6 





31 


1.7 


HD19 


2.5 


2.9 


4x900 


EG 


123.8 


2000-03-01 


KI 


3700-8500 


10 





46 


1.2 


BD26 


1.0 


1.0 


2x60 


MM, RE, JC 


137.7 


2000-03-15 


L3 


3400-5400 


7 


33 


34 


1.6 


F34 


3.0 


2.0 


800+320 


AF, AC 








5100-10500 


11 


33 


34 


1.6 


BD26 


3.0 


2.0 


800+320 


AF, AC 


158.8 


2000-04-05 


Kip 


4338-6872 


7 


66 


66 


2.3 


HD84 


1.9 


1.5 


3x900 


AF, AB 


162.8 


2000-04-09 


Kip 


3720-8685 


13 


66 


65 


2.1 


F34 


1.5 


1.5 


4x300 


RB, MB, RW 


313.1 


2000-09-06 


L3 


3400-5400 


7 


150 


154 


1.5 


BD28 


3.0 


2.0 


1800+300 


AF, RC, WL 



Table 3 — Continued 



Day a 


UT Date 


Tel. b 


Rangc c 


Res. d 


P.A.° 


Opt. P.A. f 


Air.s 


Flux Std. h 


See. 1 


Slit 


Exp. 


Observer(s) J 




(yyyy-mm-dd) 




(A) 


(A) 


(deg) 


(deg) 






(arcsec) 


(arcsec) 


(<0 










5100-10450 


n 


150 


154 


1.5 


BD17 


3.0 


2.0 


1800+300 


AF, RC, WL 


333.1 


2000-09-26 


L3 


3400-5400 


7 


180 


168 


1.5 


BD28 


2.0 


2.0 


2700 


AF, RC 








5100-7800 


7 


180 


168 


1.5 


BD17 


2.0 


2.0 


2700 


AF, RC 


516.7 


2001-03-29 


Kip 


4340-6870 


7 


60 


63 


1.7 


HD84 


1.0 


1.0 


700 


AF, AB 



a Days since discovery, 1999-10-29 (HJD 2,451,480.94), rounded to the nearest tenth of a day. 

b Ll = Lick 1-m/Nickel Reflector + spectrograph (Stone & Shields 1990); L3(p) = Lick 3-m/Kast Double Spectrograph (Miller & Stone 1993; 'p' denotes 
polarimctcr attached; 'h' denotes Hamilton echelle spectrograph used instead of the Kast); K I(p) = Keck-I 10-m/Low Resolution Imaging Spectrometer 
(LRIS [Okc ct al. 1995]; 'p' denotes polarimeter attached [Cohen 1996]); K II = Keck-II 10-m/LRIS; P = Palomar 5 m Hale/Double Spectrograph (Oke & 
Gunn 1983). 

c Observed wavelength range of spectrum. In some cases, the ends are very noisy, and are not shown in the figures. 
d Approximate spectral resolution derived from night-sky lines. 
°Position angle of the spectrograph slit. 

f Optimal parallactic angle (Filippcnko 1982) at the midpoint of the exposures. 
g Average airmass of observations. 

h The standard stars are as follows: F34 = Feige 34, BD28 = BD+28°4211— Stone (1977), Massey & Gronwall (1990); G191 = G191B2B— Massey & 
Gronwall (1990); HD19 = HD 19445, HD84 = HD 84937, BD26 = BD+26°2606, BD17 = BD+17° 4711— Okc & Gunn (1983). 

'Seeing is estimated from the FWHM of point sources on the CCD chip. 

•iAB = Aaron Barth; AC = Alison Coil; AF = Alex Filippenko; AH = Alan Harris; AQ = Andreas Quirrcnbach; AR = Adam Riess; AS = Adam 
Stanford; CG = Carol Grady; DA = David Ardila; DL = Douglas Leonard; DS = Daniel Stern; EG = Elinor Gates; GM = George Miley; JC = James 
Colbert; MB = Michael Brotherton; MM = Matt Malkan; MR = Michiel Reuland; PG = Peter Garnavich; RB = Robert Becker; RC = Ryan Chornock; 
RE = Rick Edelson; RW = Richard White; SB = Schcltc Bus; SF = Sabine Frink; WB = Wil van Breugel; WL = Weidong Li; WV = Willem de Vries. 
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As discussed in § 1, one area of uncertainty in the application of EPM is estimating the 
photospheric velocity. To test the consistency of velocities measured from different lines, we shall 
derive velocities from as many features as possible in the optical spectrum of SN 1999em during 
the photospheric phase. Our goal is to identify several unblended, weak features and compare the 
velocity derived using them with that obtained from the Fe II AA4924, 5018, 5169 lines, which have 
often been used in previous EPM studies. 

We first identified 36 distinct absorption features in the optical photospheric spectrum of 
SN 1999em, focusing mainly on the region between 4000 and 7000 A since that range was best 
covered by our spectra. We consulted previously published line identifications (IDs) as a starting 
point for major features, including Jeffery & Branch (1990), Williams (1987), Elias et al. (1988), 
McGregor (1988), Meikle et al. (1993), and Fassia et al. (1998). We then considered the theoretical 
models by Hatano et al. (1999), in which the spectral signatures of the major ions predicted to 
occur in SN spectra of all types are calculated (the results of the calculations are displayed in figures 
retrievable from the authors' web-site 20 ). To help estimate the expected relative strengths of the 
different contributing species, transition wavelengths and oscillator strengths were taken from the 
Kurucz (1996) spectral line database. We also studied the synthetic SN spectrum of progenitor 
model sl5.44.2 of Weaver & Woosley (1993) 21 , in which the non-LTE radiative transport code 
EDDINGTON (Eastman & Pinto 1993) is used to model the explosion of a 15 M Q red supergiant 
(RSG) progenitor of solar metallicity 44 days after shock breakout with an effective temperature of 
5000 K and a photosphere located at 5000 km s -1 . The optical portion of the synthetic spectrum 
is shown in Figure 8(a), with individual species responsible for many of the line features shown 
in Figure 8(b, c, d). Immediately apparent is that Fe II features dominate the spectrum below 
5500 A and consequently blend with many of the features produced by other ions. Such blended 
features are in general unreliable for simple velocity measurements since the relative strengths of the 
species making up the blend vary with photosphere location and temperature. The resulting IDs 
of the line features are shown in Figures 9 and 10, and all the likely contributors' rest wavelengths, 
configurations, and energies are listed in Table 4. The IDs for major emission-line features in the 
nebular stage are shown in Figure 11. 



20 http://www. nhn.ou.edu/~baron/papers. html. 

21 See Eastman et al. (1994) for details of the physics involved in the model's calculation; here we are primarily 
interested in its use to help identify line features. 
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I ' i ' ' I ' ' i ' I i ' ' i I ' i ' ' I ' ' i ' I i ' ' i I ' i ' 

4000 4500 5000 5500 6000 6500 7000 
Rest Wavelength (A) 

Fig. 6. — Optical flux spectra of SN 1999em spanning the first 96 days after discovery. In this 
and all figures a recession velocity of 800 km s _1 has been removed from the observed spectrum 
(§ 2.2.3). Note that the continuum shape of the spectra on days 5 and 6 has been manually adjusted 
redward of 5800 A due to irregularities introduced by the spectrograph. 
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4000 5000 6000 7000 8000 9000 10000 
Rest Wavelength (A) 



Fig. 7. — The optical and near-infrared spectral development of SN 1999em during its first 517 days 
after discovery (1999 October 29). Six plateau-phase spectra shown in Figure 6 are reproduced here 
as well, but with their full spectral range displayed. 
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Fig. 8. — (a) Synthetic spectrum of a 15M RSG progenitor with solar metallicity 44 days after 
explosion with T cff = 5000 K (model sl5.44.2 of Weaver & Woosley [1993]). (b,c,d) The relative 
strengths of line features due to Fe II, Sc II, and Ba II, respectively, in the same model. 
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Fig. 9. — Early-time flux spectrum of SN 1999em with prominent absorption features identified; 
the feature numbers corresponding to the ions listed in Table 4 are given in parenthesis. At this 
early stage, the only strong lines are due to hydrogen Balmer and He I. Weak Fe II A5169 and Si II 
A6355 features are also seen. 
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Table 4. SN II-P Photospheric-Phase Absorption Line Identifications 



Ion 


Wavelength (A) a 


log (gj) 


Feature ID b 


Designation 


E (eV) c 


HI 


3889.05 


-1.19 


1 


d 


10.20 


HI 


3970.07 


-0.99 


2 


d 


10.20 


SrII 


4077.71 


0.17 


3 


5s 2 S - hp 2 P° 


0.00 


HI 


4101.73 


-0.75 


4,3 


d 


10.20 


Fell 


4173.46 


-2.18 


5 


3d 6 ( 3 P2)4s 4 P - 3d 6 ( 5 D)4p 4 D° 


2.58 


Fell 


4178.86 


-2.48 


5 


3d 6 ( 3 P2)4s 4 P - 3d 6 ( 5 D)4p 4 F° 
3d 6 ( 1 D2)4s 2 D - 3d 6 ( 3 P2)4p 2 D° 


2.58 


Fe II 


4180.98 


-1.84 


5 


4.74 


Fe II 


4233.17 


-2.00 


6 


3d 6 ( 3 P2)4s 4 P - 3d 6 ( 5 D)4p 4 D° 
3d4s X D - 3dAp 1 D° 


2.58 


Sc II 


4246.82 


0.32 


6 


0.32 


Fe II 


4296.57 


-3.01 


7 


3d 6 ( 3 P2)4s 4 P - 3d 6 ( 5 D)4p 4 F° 


2.70 


Fell 


4303.18 


-2.49 


7 


3d 6 ( 3 P2)4s 4 P - 3d 6 ( 5 D)4p 4 D° 
hp 2 P - 6s 2 S 


2.70 


SrII 


4305.44 


-0.14 


7 


3.04 


Sc II 


4314.08 


-0.10 


7 


3d 2 3 F - 3d4p 3 D° 


0.62 


Sc II 


4320.73 


-0.26 


7 


3d 2 3 F - 3d4p 3 D° 


0.61 


Sc II 


4325.00 


-0.44 


7 


3d 2 3 F - 3dAp 3 D° 


0.60 


HI 


4340.46 


-0.45 


8,7 


d 


10.20 


Fe II 


4351.77 


-2.10 


9 


3d 6 ( 3 P2)4s 4 P - 3d 6 ( 5 D)4p 4 D° 
3d 2 3 F - 3d4p 3 F° 


2.70 


Sc II 


4374.46 


-0.44 


9 


0.62 


Fe II 


4385.39 


-2.57 


9 


3d 6 ( 3 P2)4s 4 P - 3d 6 ( 5 D)4p 4 D° 
3d 2 3 F - 3d4p 3 F° 


2.78 


Sc II 


4400.39 


-0.51 


9 


0.61 


Sc II 


4415.56 


-0.64 


9 


3d 2 3 F - 3d4p 3 F° 


0.60 


TTo TT 
re 11 


A ^ns 9Q 


9 91 
— Z.Z1 


1 1 


oil \ rZJ4S r — oil ^ Dj^ip \j 


Z.oO 


Fe II 


4522.63 


-2.03 


11 


3d 6 ( 3 F2)4s 4 F - 3d 6 ( 5 D)4p 4 D° 


2.84 


Fe II 


4549.47 


-1.75 


11 


3d 6 ( 3 F2)4s 4 F - 3d 6 ( 5 D)4p 4 D° 
3d 3 2 H-3d 2 ( 3 F)4p 2 G° 


2.83 


Ti II 


4549.62 


-0.45 


11 


1.58 


Ball 


4554.03 


0.17 


11 


6s 2 S - 6p 2 P° 


0.00 


Fe II 


4555.89 


-2.29 


11 


3d 6 ( 3 F2)4s 4 F - 3d 6 ( 5 D)4p 4 F° 


2.83 


Fe II 


4583.84 


-2.02 


11 


3d 6 ( 3 F2)4s 4 F - 3d 6 ( 5 D)4p 4 D° 


2.81 


Fe II 


4629.34 


-2.37 


12 


3d 6 ( 3 F2)4s 4 F - 3d 6 ( 5 D)4p 4 F° 
3d 21 B- 3d4p X F° 


2.81 


Sc II 


4670.41 


-0.37 


13 


1.36 


HI 


4861.32 


-0.02 


16,14 


d 


10.20 


Fe II 


4923.93 


-1.32 


17 


3d 5 As 2 6 S - 3d 6 ( 5 D)4p 6 P° 
6s 2 S - 6p 2 P° 


2.89 


Ball 


4934.08 


-0.15 


17 


0.00 


Fe II 


5018.44 


-1.22 


18 


3dHs 2 6 S-3d 6 ( 5 D)4p 6 P° 


2.89 
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Table 4. SN II-P Photospheric-Phase Absorption Line Identifications - Continued 



Ion 


Wavelength (A) a 


log (gf) 


Feature ID b 


Designation 


E (eV) c 


Sc II 


5031.02 


-0.26 


18 


3cr X D - 


- 3d4p 4 P° 


1.36 


Fe II 


5169.03 


-0.87 


19 


3d b 4s 2 b S - 


3d 6 ( 5 D)4p 6 P° 
- 3(i b ( b D)4p 4 F° 


2.89 


Fe II 


5197.58 


-2.10 


20 


3d b ( 6 G)4s 4 G 


3.23 


Fell 


5234.62 


-2.05 


20 


3d b ( 3 G)4s 4 G 


- 3d b ( 5 D)4p 4 F° 


3.22 


Sc II 


5239.81 


-0.77 


20 


- 3dAp 4 P° 


1.46 


Fell 


5276.00 


-1.94 


21 


3d b ( d G)4s 4 G 


i 1 -ik-\ / K -f— v \ a A T — 7 D 

- 3<r( b D)4p 4 F° 

- 3(i b ( 5 D)4p 4 F° 


3.20 


Fe II 


5316.62 


-1.85 


22 


3ci b ( 3 G)4s 4 G 


3.15 


Fe II 


5325.55 


-2.60 


22 


3d 6 ( 3 G)4s 4 G 


- 3<r( b D)4p 4 F° 


3.22 


Fe II 


5362.87 


-2.74 


23 


3<i b ( 3 G)4s 4 G 


- 3d b ( b D)4p 4 D° 


3.20 


Fe II 


5414.07 


-3.79 


24 


3<i b ( 3 G)4s 4 G 


- 3<i b ( 5 D)4p 4 D° 


3.22 


Fe II 


5425.26 


-3.36 


24 


3d b ( 3 G)4s 4 G 


- 3d b ( 5 D)4p 4 F° 


3.20 


Fe II 


5432.97 


-3.63 


24 


O ifi /^tt\ a 9tt 

3dr( i Yi)As Z H 
3d 2 l G - 


- 3d°(°D)4p 4 F° 

- 3dAp 4 F° 


3.27 


Sc II 


5526.79 


0.13 


25 


1.77 


Fe II 


5534.85 


-2.93 


25 


3d 6 ( 3 H)4s 2 H 


-3d 6 ( 5 D)4p 4 F° 


3.24 


Sc II 


5641.00 


-1.04 


26 


3d 2 3 P - 


- 3c%> 3 P° 


1.50 


Sc II 


5657.90 


-0.50 


26 


3d 2 3 P - 


- 3d4p 3 P° 


1.51 


Sc II 


5658.36 


-1.17 


26 


3d 2 6 P - 


- 3dAp d P° 


1.50 


Sc II 


5667.15 


-1.24 


26 


3d 2 3 P - 


- 3c%> 3 P° 


1.50 


Sc II 


5669.04 


-1.12 


26 


3<r J P - 


- 3c%> "^P 


1.50 


Sc II 


5684.20 


-1.05 


26 


3d 2 3 P - 


- 3d4p 3 P° 


1.51 


Ball 


5853.67 


-1.00 


29 


5d 2 D 


- 6p 2 P 


0.60 


ne i 


Do / O.OO 


n 7 a 


98 9Q 


ls2p 3 P° 
3s 2 S 


- Is3d 3 D 


9fl Q7 

zu.y / 


Nal 


5889.95 


0.12 


29,27 


- 3p 2 P° 


0.00 


Nal 


5895.92 


-0.18 


29,27 


3s 2 S 


- 3p 2 P° 


0.00 


Ba II 


6141.71 


-0.08 


30 


5d 2 D 


- 6p 2 P° 


0.70 


Fe II 


6147.74 


-2.72 


30 


3d 6 ( 3 D)4s 4 D 


-3d 6 ( 5 D)4p 4 P° 


3.89 


Fe II 


6149.26 


-2.72 


30 


3d 6 ( 3 D)4s 4 D 


- 3d 6 ( 5 D)4p 4 P° 


3.89 


Fe II 


6238.39 


-2.63 


31 


3d 6 ( 3 D)4s 4 D 
3d 2 3 P - 


- 3d 6 ( 5 D)4p 4 P° 

- 3dAp 3 D° 


3.89 


Sc II 


6245.64 


-0.98 


31 


1.51 


Fe II 


6247.56 


-2.33 


31 


3d 6 ( 3 D)4s 4 D 
3d 2 3 P - 


-3d 6 ( 5 D)4p 4 P° 
- 3dAp 3 D° 


3.89 


Sc II 


6279.75 


-1.21 


31 


1.50 


Sill 


6347.11 


0.30 


32 


3s 2 ( 1 S)4s 2 S 


-3s 2 ( 1 S)4p 2 P° 


8.12 


Sill 


6371.37 


-0.00 


32 


3s 2 ( 1 S)4s 2 S 
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Table 4. SN II-P Photospheric-Phase Absorption Line Identifications - Continued 



Ion 


Wavelength (A) a 


log (gf) 


Feature ID b 


Desig 


oration 


E (eV) c 


Call 


8498.02 


-1.31 


35 


3d 2 D 


- Ap 2 P° 


1.69 


Call 


8662.14 


-0.62 


36 


3d 2 D 


- Ap 2 P° 


1.69 



i 



Note. - Lines producing or contributing to absorption features in the optical spectrum of 
SN 1999em during the photospheric phase. Rest wavelengths, oscillators strengths, configurations, 
terms, and energies taken from the Kurucz (1996) spectral line database. 

a Rest wavelength (air). 

b Feature number(s) corresponding to the identifications shown in Figures 9 and 10. 
c Energy of the lower energy state of the transition. 

d For H I, the contributing closely-spaced configurations and terms are: 2p 2 P° — nd 2 D, 2s 2S — 
np 2 P°, and 2p 2 P° — ns 2 S, where n is the principle quantum number of the upper state of the 
transition (i.e., 3 for Ha, 4 for H/3, etc.). 

e For brevity, we do not list the 6 closely-spaced components of the He I A5876 multiplet, but rather 
list the mean wavelength and total oscillator strength of the line. 
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2.2.2. Line Features Associated with P-Cygni Profiles 

In all, there are five features (labeled 10, 14, 15, 27, and 33 in Figure 10) for which no convincing 
line ID is found near the rest wavelength inferred by assuming a velocity similar to that derived 
from the more readily identified lines. All five features appear at multiple epochs and follow the 
general redward drift with time seen in the other lines. We propose that three of these features (14, 
27, and 33) are actually associated with the strong P-Cygni absorptions of H/3 A4861, Na I D A5892, 
and Ha A6563, respectively, and represent the continued evolution of the low-contrast features first 
identified in the spectrum from day 7 (these features can be seen as small "notches" just blueward 
of the H/3 and Hj absorption in Figure 9; due to contamination from telluric features, the "notch" 
is not positively confirmed in Ha on day 7 - see L01). From day 30 until the end of the plateau 
these features mirror the other lines' decreasing velocity (about 2,000 km s _1 during this period of 
time), but with a velocity about 8000 km s~ x higher. It is possible that the other two unidentified 
features (10 and 15) are also high-velocity absorptions associated with strong lines: feature 10 with 
the large metal blend of Ba II, Ti II, and Fe II near A4554 (the candidate producing absorption 
most consistent with the other high- velocity features is Fe II A4584), and feature 15 with Ba II 
A4934. These identifications are less secure, however, given the relative weakness of the responsible 
lines. We note that Fassia et al. (1998) identify feature 33 with Sc II A6378. While this ID does 
yield a velocity similar to other photospheric lines, the feature is evident in neither the simulated 
spectra of Hatano et al. (1999) nor the synthetic spectrum of Weaver &: Woosley (1993) that we 
analyzed. The Kurucz (1996) line list does confirm a relatively strong Sc II feature at A6378.7, but 
it lies 7.4 eV above the ground state. For it to be visible, many other features of similar strength 
and excitation should be observed as well, and they are not. We thus favor the identification of 
this feature with Ha. 

Figure 12 shows the development of the Ha profile from day 40 until our latest observation 
on day 517. In addition to the high-velocity absorption feature already discussed (labeled A in 
Figure 12), a curious "absorption line," labeled B, develops between day 49 and day 76 and disap- 
pears by day 124. If associated with Ha, the corresponding velocity is quite similar to the velocity 
derived from weaker lines, such as Fe II A5018. A similar feature was also observed in the Type 
II-P SN 1999gi (Leonard k Filippenko 2001), as well as in SN 1988A (Turatto et al. 1993), and 
perhaps SN 1985L (Filippenko & Sargent 1986). A likely explanation is that this feature is pro- 
duced by Chugai's (1991; see also Utrobin, Chugai, & Andronova 1995) proposed mechanism for the 
anomalous blueshifted "absorption" seen in SN 1987A (i.e., the "Bochum event"; see also Phillips 
& Heathcote 1989). Under this model, the feature is actually produced by a lack of absorption by 
gas at velocities surrounding the "dip," caused by ineffective screening of continuum photons due 
to excitation stratification in the envelope. 

Finally, there is the development of a distinct "notch" in the Ha, H/3, and [Ca II] AA7291, 7324 
profiles near their rest wavelengths, which develops immediately after the end of the plateau (Fig- 
ure 12). As discussed by L01, this feature might originate within the expanding ejecta itself, 
perhaps due to an asymmetry in the line-emitting region. It is possible that it is produced by the 
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Fig. 11. — SN 1999em in the nebular phase, 333 days after discovery, with prominent emission lines 
identified; redward of A7780 the displayed spectrum is from day 313. For comparison of features, 
the spectrum from day 517 is also shown, having been scaled by a factor of 12.6 and then offset by 
-0.95; the zero-flux level is indicated (dotted line). 



- 35 - 




-10000 10000 

Velocity (km s" 1 ) 



Fig. 12. — The development of the Ha profile of SN 1999em, with days since discovery (1999-11-29 
UT) indicated. Of particular interest are the "notch" near zero velocity that is first seen on day 
124 (perhaps a hint on day 95), a "dip" (B) on the blue side of the profile on days 76 and 95, and 
a moderately strong absorption feature (A) at high velocity. The notch near zero velocity is seen 
in other strong emission lines (see L01), and the high- velocity absorption feature (A) is seen in the 
H/3 and Na I A5892 profiles as well. The plateau phase ended about 95 days after discovery, the 
same period during which features A and B disappear and the notch at zero velocity is first seen. 
Note that the spectrum from day 76 has very low resolution. 
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ejecta interacting with circumstellar material, although the lack of any other evidence suggesting 
significant circumstellar interaction (e.g., Pooley et al., 2001) makes us doubt this possibility. The 
notch seems to diminish in strength with time, and it may not be present on day 517, although 
the low signal-to-noise ratio (S/N) of that spectrum renders it insensitive to low-contrast features. 
It is also interesting to note the many other features that complicate the Ha absorption profile as 
SN 1999em makes the transition out of the plateau epoch, particularly evident in the day 95 spec- 
trum; these features are also seen in a spectrum of SN 1999gi at a similar epoch; in fact, the exact 
correspondence of all of the line featurs between these two SNe II-P is quite remarkable (Leonard 
& Filippenko 2001). 



2.2.3. Inferring Photospheric Velocity from Absorption Lines 

From the lines identified in Figure 10, we find four weak, unblended spectral features that 
should serve as good photospheric velocity indicators to compare with the velocity determined by 
the stronger Fe II AA4924, 5018, 5169 lines: Fe II A4629 (feature 12), Sc II A4670 (feature 13), Fe II 
A5276 (feature 21), and Fe II A5318 (feature 22). To translate the measured wavelength minima 
into velocities we must remove the redshift of SN 1999em and assign rest wavelengths to all of the 
line features. The NASA/IPAC Extragalactic Database (NED) recession velocity of NGC 1637 is 
717 km s _1 . However, since the host galaxy is somewhat inclined (i = 32°) and the location of the 
SN within the spiral arm is not known, the NED velocity may differ significantly from the actual 
velocity of the SN. One way sometimes used to infer the recession velocity of a SN is to measure 
the velocity of bright H II regions on both sides of the SN along the slit, and interpolate to the 
position of the SN. When this is done for SN 1999em on the frames that contain bright H II regions, 
however, values from 690 to 900 km s _1 result; in other words, we find the inferred velocity to be 
a function of which H II regions fall along the slit, which is determined by the slit's position angle. 
As shown in Figure 14, the Na I D interstellar absorption from NGC 1637 occurs over a range from 
about 700 to 820 km s _1 . We are thus left with some uncertainty in the exact recession velocity 
of SN 1999em. This uncertainty ultimately leads to a small systematic uncertainty in our derived 
EPM distance to SN 1999em. For our analysis, we assign a recession velocity to SN 1999em of 800 
km s _1 and a la uncertainty of 50 km s _1 . 

Assigning a single rest wavelength to each of the features identified in Figures 9 and 10 is 
complicated by the blending that exists in many of the lines, and simply taking the rest wavelength 
of the dominant line may lead to inaccuracies. When a feature consists of a closely spaced multiplet 
of a single ion, we estimate the rest wavelength as 

x _ E"=i A * x (jjfh n <v, 
Ei=i(5/)i 



where i is the line number and (gf) is the oscillator strength. Since it has been previously used 
as a photospheric velocity indicator (e.g., Schmidt et al. 1994a), we note that the 6 lines of the 
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multiplet of Sc II comprising feature 26 yield A = 5665.6 A, which is significantly different from the 
rest wavelength of the strongest component, A = 5657.9 A; using A = 5657.9 A, as was evidently 
done in previous studies, would therefore result in derived velocities that are ~ 400 km s -1 lower 
than is appropriate for this feature. For features that contain blends of lines from different ions (or 
different energy levels within the same ion), the assigned rest wavelength was guided by estimates 
of the relative strengths of the individual lines during the photospheric phase. In many cases, 
though (i.e., Fe II A4924 and Sc II A5527), the relative strengths of the individual lines likely vary 
considerably, diminishing their reliability as simple photospheric velocity indicators. 

To assign an uncertainty to the inferred velocities, we attempt to quantify the error resulting 
from uncertainty in both the wavelength scale and the measurement process itself. To gauge the 
uncertainty in our wavelength scale, we cross-correlated the region around the strong Na I D 
AA5890, 5896 host galaxy interstellar absorption line in every low-resolution spectrum with this 
region in our day 95 spectrum, which has an excellent S/N. The la spread of wavelength offsets 
was found to be 0.6 A, which we adopt as the uncertainty in our derived wavelength scale for every 
spectrum. 

To measure the wavelength of maximum absorption and derive its measurement uncertainty, 
we first estimated the local "continuum" shape by fitting a straight line to points near the two local 
flux peaks on either side of the absorption feature. We then normalized the absorption feature by 
dividing it by the interpolated continuum fit. To estimate the wavelength of minimum flux in the 
normalized line profile, we then fit a series of third-order Savitsky-Golay smoothing polynomials 
(Press et al. 1992) to the data with widths ranging from 1/2 to 3/2 the FWHM of the line profile 
and then averaged the fits; this was found to produce more robust results than using a single 
smoothing function alone. The minimum of the smoothed spectrum was then found, interpolated 
to the nearest 0.01 A. To estimate the measurement error, we added random noise at the level 
seen in the data to the smoothed spectrum and used an automatic routine to find the absorption 
minimum in 1000 sets of noisy artificial data. The standard deviation of these measurements added 
in quadrature to the 0.6 A uncertainty in our wavelength scale is the final uncertainty reported 
for the line's absorption minimum. We do not attempt to quantify the systematic uncertainty in 
a feature's assigned rest wavelength, but caution that for some blends this systematic error may 
dominate the total (statistical and systematic) uncertainty. 

We characterize each feature's strength by measuring its "depth," defined as 

d=(fc~ /min)//c, (13) 

where f m i n is the flux at the line's minimum and f c is the value of the interpolated continuum flux at 
the location of the line's minimum. The measured depth in some cases depends significantly on the 
strength of adjacent features since this affects the location of the points used to fit the normalizing 
continuum. In general, though, d should be a reasonably good indicator of line strength. The rest 
wavelengths adopted for all features and the measured velocities, velocity uncertainties, and line 
depths are given in Table 5. 



Table 5. SN 1999em: Absorption Line Measurements 
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40) 

40) 
(40 

40) 

40) 

50) 

30) 
(40 

(50) 



3370 (30) 
3100 (30) 
3470 (30) 

11640 (30) 
4830 (80) 



0.26 
0.12 
0.62 
0.17 
0.14 
0.13 
0.07 
0.32 
0.40 
0.59 
0.59 
0.09 
0.13 
0.19 
0.04 
0.11 
0.32 
0.35 
0.02 

0.59 
0.19 
0.25 

0.15 
0.73 



3290 (330) 0.27 



2790 (160) 
2290 (70) 
2310 60) 
10550 (170 
[4765.0 (1.6 
2310 (70) 
2550 (60) 
2790 (50) 
2900 (130) 



2240 (110) 
2460 (70' 
2580 (60 



3010 (90) 
2610 (70) 
2800 (50) 

11300 (90) 
5440 (120 
4150 (260 
3660 (160 



0.69 
0.22 
0.21 
0.12 
0.25 
0.28 
0.50 
0.70 
0.70 



2620 (290) 0.42 



0.20 
0.44 
0.52 



2820 (250) 
[4422.6 (1.7)] 



2010 (90) 
1850 90) 
11240 (100 
[4755.5 (1.6 
2060 (50) 
2240 (70) 
2500 (60) 
2660 (110) 
2100 (110) 
2110 (80) 
1950 (40 
1570 (50 
1850 (40 
2140 (50 
2360 (30 



0.34 
0.31 
0.11 
0.05 
0.39 
0.61 
0.76 
0.64 
0.16 
0.25 
0.35 
0.14 
0.24 
0.51 
0.61 



0.33 
0.22 



2450 ( 


500) 


0.66 


2350 (80) 


0.66 


1900 ( 


100) 


0.38 








2060 


(80) 


0.23 








10110 


(440 


0.10 


10700 


90 


0.13 


[4775.2 (1.9)] 


0.19 


[4783.2 


0.9) 


0.12 


2120 (160) 


0.34 


2010 ( 


40) 


0.42 


2330 




0.61 


1830 ( 


50) 


0.52 


2600 




0.73 


2030 ( 


40 


0.58 


2660 ( 


190 


0.64 


1960 ( 


40) 


0.62 


2110 ( 


200 


0.17 








2120 (160) 


0.21 








2030 (200) 


0.32 


1710 (30) 


0.36 


1970 


(80) 


0.22 


1340 ( 


40) 


0.42 


2180 


60) 


0.47 


1450 ( 


40) 


0.17 


2250 


(60) 


0.58 


1790 ( 


40) 


0.40 



0.70 


3240 (30) 


0.76 


3140 (60) 


0.76 


5620 (40) 


0.89 


0.41 


2640 (40) 


0.52 


2410 (80) 


0.53 


2720 (30 


0.60 


0.38 


2350 (40) 


0.46 


2190 (50) 


0.45 


1770 (40) 


0.31 


0.17 














0.63 


5100 (610) 


0.65 


5600 (150) 


0.64 


4440 (50) 


0.83 


0.85 






2890 (130) 


0.82 



to 



0.54 



Table 5 - Continued. 





Day 138 


Day 159 


Day 163 


Day 313 


Day 333 


Line 


Vel. ' Depth 


Vel. Depth 


Vel. Depth 


Vel. Depth 


Vel. ' Depth 



1 

2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 



1800 



(140) 0.62 1950 



2190 (60) 
1790 (100) 
1910 (130 
1430 (220) 



1510 (300) 
1260 (80' 
1630 (70 



5330 (50 
2660 (60 
1660 (170) 



0.75 
0.53 
0.54 
0.67 



1500 (130) 0.41 



0.40 
0.18 
0.44 



0.91 
0.58 
0.30 



(190) 



2360 (70 
1510 (40 
1780 (40) 
960 (40) 



5650 (160) 
2730 (60) 
1610 (30) 



0.61 



1620 



(140) 



0.85 
0.53 
0.59 
0.67 



10500 (100) 0.06 



0.92 
0.49 
0.39 



1150 (70) 
1290 (60' 
1500 (90 



5480 (130) 
2770 (80) 
1630 (50) 



0.57 1290 



10760 (70) 0.12 10780 (110) 0.11 



2490 (60 
1530 (50 
1780 (50) 
1000 (120) 



1430 (80) 0.38 



0.41 
0.17 
0.39 



0.91 
0.46 
0.38 



(390) 0.67 



1800 (190) 
[4436.9 (3.7)] 
2280 (190) 



4840 (430) 
4540 (270) 
1770 (290) 



0.95 
0.54 
0.61 



4060 (460) 
4930 (150) 
1700 (80) 



0.39 
0.28 
0.50 



0.84 


2190 


240 


1.05 


2370 




0.86 


0.48 


1500 


100 


0.58 


1520 




0.36 


0.58 


1780 


;2so; 


0.66 


1910 




0.60 


0.65 


1910 


,180, 


0.30 


1910 




0.24 



2450 (170) 0.42 2260 (90) 0.48 



0.91 
0.53 
0.60 



3740 (110) 0.88 3580 (30) 0.92 3630 (90) 0.92 3460 (70) 0.98 3520 (30) 0.96 



Table 5 - Continued. 





Day 138 


Day 159 Day 163 


Day 313 


Day 333 


Line 


Vel. Depth 


Vel. Depth Vel. Depth 


Vel. Depth 


Vel. Depth 


35 
36 


2980 (40) 0.93 
2970 (40) 0.63 


2910 (180) 0.88 

2660 (50) 0.65 


2440 (190) 0.91 
2250 (110) 0.84 





Note. - Line velocities for absorption features measured directly from spectra of SN 1999em with 
estimated la statistical uncertainty shown in parentheses, shown for four example days only. Uncertainty 
includes measurement error and uncertainty in the wavelength scale, but does not include uncertainty in the 
rest wavelength estimation, which may be significant for blended features. All spectra were deredshifted by 
800 km s _1 prior to the velocity determination. Note that the emission-dominated nebular spectrum from 
day 517 is not included in the tabulation. 

a Day since discovery, 1999 October 29 (HJD 2,451,480.94), rounded to the nearest day. 

b Line number corresponding to feature indicated in Figures 9 and 10, and identified in Table 4. 

c Main ion responsible for feature. Colon indicates significant blending of multiple ions; see Table 4 for 
components of blend. Velocities derived using blended lines may not yield accurate photospheric velocities, 
since the relative strengths of the components may change with time. 

d Rest wavelength (air) adopted for line feature. In the case of blended lines, a weighted average of the 
blended components is used; see text for details. 

e Velocity measured in units of km s^ 1 , with uncertainty listed in parentheses. For line features lacking 
identification, the measured rest wavelength and uncertainty are given in brackets. 

f Depth is defined as: d = 1 — /mim where f m \ n is the normalized flux level at the location of the line's 
minimum with the continuum set to 1. It should be noted that the measured depth in some cases may 
depend on the strength of adjacent features since this affects the location of the points used to fit the 
normalizing continuum. 

g To maintain consistency with the method traditionally used to determine the photospheric velocity 
by averaging the velocities derived from Fe II AA4924, 5018, 5169 line measurements, we adopt the rest 
wavelength of Fe II A4924 to determine this feature's velocity. The consistently lower velocity derived 
compared with the other Fe II lines of the multiplet, however, suggests significant blending with Ba II 
A4934. 

h Becomes weak and probably blended with Na I D by day 10. 
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2.3. The Reddening of SN 1999em 

SN 1999em suffers from minimal Galactic reddening, E(B—V)mw = 0.04 mag (Schlegel, 
Finkbeiner, & Davis 1998). From theoretical fits to the early-time continuum shape, Baron et 
al. (2000) estimate the total (Galactic + host galaxy) reddening of SN 1999em to be 0.05 < 
E(B—V)tot ^ 0.10 mag, with an upper limit of E(B— V)tot < 0.15 mag; the most successful 
fits seem to favor E(B—V)tot = 0.10 mag. We checked the consistency of this result using two 
techniques. 

First, we compared the BVI color of SN 1999em with that of SN 1987A at the moment the 
He I A5876 line disappears from both objects' spectra. Since exciting a significant amount of He I 
to the high energy level (E = 20.97 eV; see Table 4) needed to produce the A5876 line is very 
temperature sensitive (generally greater than about T c g- = 10, 000 K; Hatano et al. 1999; Branch 
1987; see also Figure 1), we expect this line to disappear when the photosphere is at a well-defined 
temperature. Therefore, the intrinsic color of the two objects should be similar at this early phase, 
before significant line blanketing from metal lines begins. We assume the reddening of SN 1987A 
to be E{B—V) = 0.15 =L 0.05 mag (see discussion in Arnett et al. 1989) and compare its color 
at this phase with that of SN 1999em. For SN 1999em, the He I A5876 line becomes too weak to 
measure 10 ± 1 days after discovery, roughly corresponding with the epoch when the Fe II A5169 
absorption first appears (strengthening the argument that at this phase the optical spectrum is 
nearly a blackbody, with little line blanketing). At that epoch, its BVI magnitudes are 14.01, 
13.85, and 13.49, respectively, interpolated from the values given in Table 2. Phillips et al. (1988) 
estimate that the He I A5876 line disappeared in the spectrum of SN 1987A on JD 2,446,854, just 
four days after discovery. From the photometry of Hamuy et al. (1988), the BVI magnitudes for 
SN 1987A on that date are 4.83, 4.45, and 4.01, respectively. Using the reddening law of Savage & 
Mathis (1979) and minimizing the color difference between the two SNe at this epoch, we derive a 
reddening of E{B— V) = 0.03 ± 0.05 mag for SN 1999em, in good agreement with the estimate of 
Baron et al. (2000). 

A second method sometimes used to provide a reddening estimate for SNe II-P (e.g., Schmidt 
et al. 1993) is to assume that the early-time color evolution of all SNe II-P is similar, and compare 
the early color evolution of a SN II-P with that observed in SN 1969L, a Type II-P believed to 
suffer nearly zero host-galaxy reddening (Kirshner & Kwan 1974; Schmidt et al. 1992), and only 
E(B—V)mw = 0.06 mag (Schlegel et al. 1998). Using the time of explosion for SN 1999em derived 
in § 3 (5.3 ±1.4 days before discovery), Figure 13 shows that at early times, SN 1999em and SN 
1969L have a very similar color. This suggests similarly low reddenings for the two events. We do 
note, however, that the B — V color evolution of SN 1999em at later times deviates substantially 
from that of SN 1969L, and may indicate that SN 1999em's subsequent color evolution differed 
somewhat from the norm. Although limited by the small number of observations and rather poor 
temporal coverage, Schmidt et al. (1992) demonstrate that the B — V color evolution of the 4 
SNe II-P in their sample is quite similar to that of SN 1969L during the recombination phase 
(none has photometric data taken prior to 14 days after the explosion, however). It is possible 
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Fig. 13. — The temporal evolution of the B — V color of SN 1999em compared with that of the 
prototypical Type II-P event, SN 1969L, which is believed to be nearly unreddened. The time of 
explosion for both events is derived through the EPM analysis (see § 3 for SN 1999em and Schmidt 
et al. [1992] for SN 1969L). The good agreement between the early-time color evolution of the two 
events suggests that SN 1999em might also suffer little reddening. 
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Fig. 14.— Interstellar Na I D (AA5890, 5896) absorption by gas in NGC 1637 along the line of sight 
to SN 1999em. 
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that the large deviation between the later color evolution of SN 1999em and SN 1969L indicates a 
real difference between SN 1999em and this prototypical SN II-P. We shall further investigate the 
possibility of unusual photometric behavior for SN 1999em in § 4.4. 

A final indication that there may be little reddening for SN 1999em comes from the low value 
of the interstellar polarization implied by the spectropolarimetry (L01). One argument for higher 
extinction, however, comes from the strength of the Na I D (AA5890, 5896) interstellar absorption 
lines. From the echelle spectrum taken 4 days after discovery (Figure 14; a detailed analysis of the 
Na I D interstellar lines will be the subject of a future paper), we measure W\ (D2; A5890) = 0.86 
A and W\ (Dl; A5896) = 0.65 A, for a total equivalent width of W x (Na I D) = 1.51 A (cf. Jha et 
al. 1999). If we adopt the rough correlation derived by Barbon et al. (1990) between W\ (Na I D) 
and reddening, 

E{B—V) = 0.25W"A(Na I D), (14) 

we would predict E(B— V) ~ 0.38 mag for SN 1999em, significantly greater than the extinction 
implied by the other methods. The relation given by equation (14) is known to have very large 
scatter (likely due to variations in the dust-to-gas ratios among galaxies as well as the marginal 
ability of Na I D gas to serve as a tracer of the hydrogen gas column; saturation effects can also 
play a large roll; see, e.g., Munari & Zwitter 1997; Issa, MacLaren, & Wolfendale 1990), however, 
and we are inclined to favor the evidence supporting a lower total reddening; in § 4.1 we shall also 
see that the EPM distances derived for E{B— V) = 0.38 mag using the different filter combinations 
are very inconsistent with each other, further arguing against such high extinction. We therefore 
adopt the reddening preferred by Baron et al. (2000) of E{B—V) = 0.10 mag and assign it an 
uncertainty of ± 0.05 mag. The strength of the Na I D lines, however, does make us treat this 
value with some caution. 



2.4. The Photospheric Velocity and Color of SN 1999em During the Plateau 

From Figure 3 we see that the plateau phase of SN 1999em lasted until about 95 days after 
discovery, and our application of EPM will therefore be restricted to times somewhat earlier than 
the drop off the plateau. Figure 15 shows the velocity inferred from the flux minima of the first 
3 hydrogen Balmer lines and Fe II A5169 during the first 130 days after discovery. As expected 
(§ 1), the Ha and H/3 lines generally yield velocities significantly higher than the weaker H7 and 
Fe II A5169 lines. Curiously, after about day 40, the H/3 line yields velocities similar to, or even 
less than, those derived from Fe II A5169. In the day 76, 95, and 96 spectra the Ha absorption 
profile becomes corrupted by the excess emission described in § 2.2.2 (Figure 12), and the inferred 
velocity increases dramatically. 

Traditionally, the Fe II AA4924, 5018, 5169 line features have been used to estimate photospheric 
velocity during the plateau phase of SNe II-P, with the occassional use of Sc II A5527 and Sc II 
A5658 as well. At very early times, before the iron lines become available, it has been customary to 
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Fig. 15. — Line velocity determined from the P-Cygni absorption minima of 4 different features 
during the first 130 days after discovery of SN 1999em. 
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either use the higher-order hydrogen Balmer lines (i.e., H7 and H5) or else to scale the Ha and H/3 
velocities using the ratio of their velocities to the iron features when the iron features first appear. 
Here, we shall follow two approaches. First, we derive velocities following the usual methods for 
determining photospheric velocity (^strong)) relying mainly on the Fe II AA4924, 5018, 5169 lines 
when they (or, at least, one of them) become available (day 7 onwards), and the higher-order 
hydrogen Balmer lines at earlier times. Second, we estimate photospheric velocity by using the 
weakest unblended features identified during each epoch (v wea k)- We shall then compare the EPM 
distance derived using each set of velocities. 

Figure 16 shows the velocities measured for features observed during the first 11 days after 
discovery of SN 1999em. Through day 5, Ha and H/3 yield similar velocities. He I A5876 and H7 
yield similar velocities through day 6, but significantly lower than those measured by Ha and H/3. 
From day 7 onwards, the velocities inferred from the weak He I A5876 and H(5 lines are somewhat 
lower than those measured using Fe II A5169 and H7. Our estimates of ^strong and v wea k, then, 
will be the same for days 2 through 6, but will differ from day 7 onwards. For days 2 through 6, 
we take the weighted average of the He I A5876 and H7 velocities. From day 7 onwards, throng is 
always determined by the available Fe II AA4924, 5018, 5169 lines, and w wea k is derived from the 4 
weak, unblended features described in § 2.2.3 (features 12, 13, 21, and 22 in Figure 10). 

Figure 17 shows the velocities derived from these features compared with those derived from 
the weighted average of the Fe II AA4924, 5018, 5169 lines from day 20 to day 100. Clearly, these 
weaker features yield systematically lower velocities than those derived from the stronger Fe II lines, 
and the difference appears to grow with time. To determine v W eak, then, we took the weighted mean 
of the available weak lines at each spectral epoch. Table 6 lists the ^strong and f wea k values adopted 
for each epoch; the reported uncertainty is the uncertainty of the weighted mean of the lines added 
in quadrature to the la spread in the velocities of the individual lines. 

Figure 18 compares ^strong and f wca k during the plateau phase of SN 1999em. Since lower 
photospheric velocity estimates will result in a smaller derived distance (equation 3), it is useful to 
see whether such a velocity offset is ubiquitous among SNe II. Fig ure 19 shows the ratio of fweak 
to ^strong as a function of time in the spectra of 6 other SNe II that had the high S/N needed to 
make the measurements of these weak features. The good agreement of the data from other SNe II 
with SN 1999em suggests that EPM distances derived using the strong iron lines as photospheric 
velocity indicators will generally yield distances 5 — 10% larger than those using the weaker lines. 

In order to derive an EPM distance to SN 1999em, we need its BVI magnitude at each of 
the spectral epochs. The spectral observations of SN 1999em listed in Table 3 are not, in general, 
coincident with the BVI photometric observations listed in Table 2. We therefore used a linear 
interpolation from the nearest photometric observations before and after the spectral observations. 
The reported photometric uncertainty represents the quadrature sum of the uncertainty in the 
interpolation (i.e., uncertainty in the location of the interpolated point due to the uncertainty in 
the slope of the line connecting the two photometric data points) with the photometric uncertainty 
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Fig. 16. — The velocity derived from absorption lines in the early-time spectrum of SN 1999em. 
Note that from day 7 onwards the velocity found for Fe II and H7 is somewhat higher than the 
weaker RS and He I A5876 lines suggest. 
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Fig. 17. — Temporal evolution of the velocity inferred from weak, unblended features compared 
with that derived from the stronger Fe II AA4924, 5018, 5169 absorptions. The weaker features 
give consistently lower velocities. 
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Table 6. Quantities Used in the EPM Analysis of SN 1999em 





^strong 


^weak 




B(a B ) 




V(a v ) 






Day a 


km s _1 


km s _1 




mag 




mag 




mag 


1.9 


9276(201) 


9276(201) 


13 


sni 8 


0971 


13 


808^ 
oUou 


'0 09/11 ^ 


13 




'0 0987^ 


2.9 


9221(192) 


9221(192) 


13 


7QM 


09^9^1 
U.UZoZ J 


13 


7887 


'Ci 09/1^ 


13 


ODZU 


098^ 


3.9 


8581(70) 


8581(70) 


1 1 

13 


891 A 


031 1 


1 o 

13 


7Q97 


'Ci 0901 ^ 


1 o 

13 


^^79 


09fi1 


4.9 


8310(77) 


8310(77) 


1 o 

13 






1 o 

13 


7877 


'Ci 01 ^\ 


1 o 

13 


^41 1 

0411 


094^ 


5.9 


8314(171) 


8314(171) 


13 


OOUO 




13 


8DQ^ 

ouyo 


'Ci 01 7&\ 


13 


^ZLQ 


09Q9"l 


7.0 


8026(178) 


7372(110) 


13 






13 


8371 


'0 01 AC\\ 


13 




C\ 09Q0"l 


7.9 


8080(402) 


6809(372) 


13 




o ni qi ^ 


13 


83Q8 


'Ci 01 30"\ 


13 


ZLQ77 


09^1 ^ 


8.9 


7642(271) 


6759(317) 


13 


Q8^7 


o ni m ^ 


13 


8^8 

OOUO 


'Ci 0931 "\ 


13 


^079 

OKJ t Z 


091 1 ^ 


10.2 


7221(155) 




14 


uuoy 




13 


8481 

O^Ol 


'Ci CiACiA} 


13 




n OzL^^i^ 

U.U'iOtJ ) 


11.2 


TO 1 r / (\ A \ 

7215(94) 




1 A 

14 


0354 


;0.0559) 


1 o 

13 


8389 


U0428) 


1 '} 

13 


4702 


;0.0404) 


20.9 


K f* 1 1 /1 1 f \ 

5611(115) 


5315(174) 


14 


4667 


U0230) 


13 


8579 


;0.0290) 


13 


4054 


;0.0270) 


23.9 


5110(122) 


4887(403) 


14 


5791 


;0.0262) 


13 


8804 


;0.0323) 


13 


3818 


U0309) 


25.9 


4902(169) 


4872(305) 


14 


6532 


U0209) 


13 


8953 


;0.0197) 


13 


3663 


;0.0253) 


29.9 


4623(52) 


4358(134) 


14 


7882 


;0.0371) 


13 


9358 


;0.0250) 


13 


3490 


;0.0350) 


33.9 


4310(79) 


4047(86) 


14 


8673 


U0250) 


13 


9272 


;o.o2oo) 


13 


3118 


;o.o2io) 


37.1 


3809(89) 


3581(107) 


14 


9257 


U0220) 


13 


9209 


;0.0222) 


13 


2919 


;0.0275) 


39.9 


3785(86) 


3596(86) 


14 


9741 


;o.oi80) 


13 


9298 


;o.oi90) 


13 


2845 


U0320) 


43.9 


3836(69) 


3550(140) 


15 


0025 


;0.0467) 


13 


9123 


;0.0369) 


13 


2226 


;o.o5oo) 


46.9 


3747(103) 


3443(64) 


15 


0541 


;o.o3oo) 


13 


9270 


;0.0350) 


13 


2200 


;0.0390) 


48.8 


3500(141) 


3242(34) 


15 


0800 


U0344) 


13 


9387 


;0.0281) 


13 


2148 


;0.0362) 


75.9 


2714(164) 


2312(78) 


15 


4031 


U0373) 


14 


0396 


U0271) 


13 


2176 


;0.0333) 


94.8 


2430(190) 


1965(87) 


15 


7184 


U0227) 


14 


2587 


;0.0224) 


13 


3670 


U0328) 


95.8 


2504(170) 


2012(106) 


15 


7386 


;o.o2io) 


14 


2770 


;0.0220) 


13 


3795 


;0.0320) 



Note. - Although only data up through day 76 were included when deriving the 
EPM distance, we list the values of the measured parameters through the early part of 
the transition to the nebular phase for completeness. 

a Days since discovery, 1999-10-29 UT (HJD 2,451,480.94). 
b Velocity derived largely from Fe II AA4924, 5018, 5169 lines. 
c Velocity derived from weak, unblended lines; see text for details. 
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Fig. 18. — Comparison of photospheric velocity derived for SN 1999em by using the average of the 
strong Fe II AA4924, 5018, 5169 lines (^strong) or the weaker unblended features (f W eak) shown in 
Figure 17. 
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Fig. 19. — The temporal evolution of the ratio of the velocity derived from weak unblended lines to 
that derived from the strong Fe II AA4924, 5018, 5169 features for 7 SNe II. Note that the relative 
phases of the SNe are somewhat uncertain since some may have been discovered significantly after 
explosion. 
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of the nearest photometric point. TheBVI magnitude and uncertainty adopted for SN 1999em at 
the times of spectral observations are given in Table 6. 

3. The EPM Applied to SN 1999em 

Following the procedure outlined in § 1, we derive the theoretical angular size (6), flux dilution 
factor (£), and color temperature of SN 1999em for the three filter combinations BV, BVI, and VI, 
and list the results in Table 7. The uncertainty given for 6, £, and T at each epoch is the la spread 
found for each parameter derived from 1000 simulated sets of photometry characterized by the 
photometric uncertainty listed in Table 6. 

To determine the distance to SN 1999em for each filter combination we first need to decide how 
to weight the data, and then find the best-fitting line in the plot of t vs. 6/v (equation [3]). Since 
the uncertainty of 6/v increases in direct proportion to the value of 6/v, a weighting scheme based 
on the calculated uncertainty in each point is necessarily biased towards the earlier data. On the 
other hand, using the method of least squares with equal weighting is susceptible to undue influence 
by deviant points. We therefore opt for the more robust fitting method of using the criterion of 
the least absolute deviation for equally weighted data points (i.e., 'medfit' from Numerical Recipes; 
Press et al. 1992). Since SN 1999em begins to drop off the plateau near day 95 after discovery, 
we do not consider the data from day 95 onwards in the EPM analysis. The distances and times 
of explosion resulting from the 3 filter combinations and two photospheric velocity estimates are 
given in Table 8, along with uncertainties derived from 1000 simulated sets of data. Figures 20 
through 25 show the relevant plots for each of the six distance estimates. 

While the formal uncertainty of each derived distance listed in Table 8 is only a few percent, 
discrepancies among the distances derived using the dilution factors given for the various filter com- 
binations run as high as 19%. It is therefore clear that systematic, not statistical, error dominates 
the uncertainty in the EPM distance to SN 1999em. In general, Dyj > Dbvi > Dbv- We note 
that the 50 km s _1 uncertainty in the redshift of SN 1999em (§ 2.2.3) contributes an uncertainty of 
only ~ 2% to the derived distance. As expected, distances derived using the stronger iron features 
(^strong) are systematically about 10% larger than those derived using the weaker lines (v wea k)- For 
the weak lines, the (unweighted) average of the distances derived from the 3 filter combinations is 
Avcak = 8.24 ± 0.56 Mpc, with t Q = 5.30 ± 1.45 days before discovery; the reported uncertainties 
are the standard deviations of the three individual data points from the mean value. 22 Similarly, 
for the strong lines, -D s trong = 9.03 ± 0.78 Mpc with t Q = 5.69 ± 1.70 days before discovery. 

Since we believe that photospheric velocity is slightly overestimated by using the traditional 
Fe II AA4924, 5018, 5169 lines, we conclude that the most accurate distance to SN 1999em is that 



22 We note that using the method of least squares to determine the best-fitting slope and intercept of the line 
described by equation 3 yields D wca k = 8.54 ± 0.39 Mpc, with t = 6.20 ± 0.92 days before discovery. 



Table 7. Quantities Derived from the EPM Analysis of SN 1999em 





Angular Size (10 s km Mpc _1 ) a 




Dilution Factor 




Temperature (K) 


Day b 


Obv Obvi Ovi 


Cbv 


Cbvi 


Cvi 


Tbv Tbvi Tvi 



I. 9 
2.9 
3.9 
4.9 
5.9 
7.0 
7.9 
8.9 
10.2 

II. 2 
20.9 
23.9 
25.9 
29.9 
33.9 
37.1 
39.9 
43.9 
46.8 
48.8 
75.9 
94.8 
95.8 



6.13(1.01) 
6.54(0.94) 
7.14(1.09) 
7.91(0.81) 
8.25(0.97) 
8.55(0.83) 
9.18(0.67) 
9.59(0.82) 
10.52(1.70) 
11.45(1.79) 
16.46(0.38) 
16.64(0.36) 
16.73(0.23) 
16.75(0.32) 
17.16(0.27) 
17.51(0.30) 
17.65(0.28) 
18.06(0.61) 
18.18(0.59) 
18.18(0.52) 
19.47(0.79) 
18.94(0.67) 
18.82(0.65) 



7.57(0.50) 
7.79(0.50) 
8.12(0.51) 
8.56(0.47) 
8.99(0.50) 
9.45(0.48) 
10.12(0.37) 
10.31(0.31) 
10.81(0.73) 
11.28(0.79) 
15.00(0.31) 
15.69(0.34) 
16.06(0.31) 
16.52(0.29) 
17.00(0.20) 
17.29(0.24) 
17.40(0.24) 
17.94(0.36) 
18.00(0.32) 
18.04(0.29) 
18.14(0.23) 
16.94(0.22) 
16.84(0.22) 



8.46(0.70) 
8.59(0.69) 
8.76(0.66) 
9.01(0.54) 
9.50(0.68) 
10.07(0.67) 
10.78(0.56) 
10.85(0.58) 
11.12(1.13) 
11.34(1.09) 
13.00(0.76) 
13.82(0.76) 
14.34(0.59) 
15.22(0.77) 
15.84(0.47) 
16.15(0.58) 
16.40(0.65) 
17.39(1.03) 
17.60(0.85) 
17.82(0.76) 
18.72(0.62) 
17.99(0.51) 
17.93(0.51) 



0.411( 
0.404( 
0.394( 
0.384( 
0.380( 
0.376( 
0.37K 
0.369( 
0.367( 
0.367( 
0.598( 
0.699( 
0.776( 
0.914( 
1.06K 
1.180( 
1.258( 
1.352( 
1.43K 
1.462( 
2.002( 
2.259( 
2.264( 



0.020 
0.017 
0.016 
0.009 
0.010 
0.007 
0.004 
0.003 
0.007 
0.010 
0.040 
0.050 
0.040 
0.071 
0.056 
0.058 
0.052 
0.123 
0.104 
0.101 
0.120 
0.086 
0.083 



0.429(0.005) 
0.427(0.005) 
0.424(0.005) 
0.420(0.004) 
0.417(0.003) 
0.414(0.002) 
0.412(0.001) 
0.412(0.000) 
0.412(0.001) 
0.412(0.002) 
0.467(0.008) 
0.497(0.011) 
0.520(0.011) 
0.565(0.015) 
0.603(0.011) 
0.632(0.014) 
0.652(0.014) 
0.690(0.025) 
0.714(0.023) 
0.728(0.020) 
0.890(0.030) 
0.988(0.025) 
0.994(0.022) 



0.444(0.005) 
0.444(0.005) 
0.443(0.004) 
0.441(0.003) 
0.439(0.003) 
0.437(0.002) 
0.436(0.001) 
0.436(0.001) 
0.436(0.002) 
0.436(0.002) 
0.440(0.004) 
0.445(0.006) 
0.450(0.005) 
0.460(0.009) 
0.466(0.007) 
0.470(0.008) 
0.474(0.009) 
0.486(0.018) 
0.491(0.016) 
0.496(0.015) 
0.531(0.017) 
0.560(0.018) 
0.562(0.018) 



17488(2121) 
16809(1726) 
15780(1739) 
14709(1063) 
14177(1192) 
13703(945) 
13012(667) 
12529(716) 
11723(1481) 
10990(1465) 
6523(249) 
6005(215) 
5706(141) 
5286(184) 
4948(114) 
4725(100) 
4598(81) 
4460(172) 
4354(132) 
4315(125) 
3789(93) 
3605(56) 
3602(54) 



14050(646) 
13834(625) 
13384(637) 
12882(524) 
12361(480) 
11829(416) 
11250(282) 
11019(229) 
10640(530) 
10294(605) 
7612(151) 
7110(162) 
6822(127) 
6378(119) 
6085(77) 
5899(84) 
5778(79) 
5582(120) 
5472(100) 
5410(86) 
4857(82) 
4609(57) 
4596(49) 



12150(798) 
12094(773) 
11915(701) 
11697(531) 
11185(614) 
10641(527) 
10135(394) 
10032(431) 
9886(823^ 
9779(781)ot 
8820(439)7 
8375(403) 
8108(278) 
7649(327) 
7442(201) 
7348(235) 
7238(248) 
6957(378) 
6854(319) 
6756(271) 
6242(207) 
5925(169) 
5900(167) 



Note. — Although only data up through day 76 were included when deriving the EPM distance, we list the values of the measured 
parameters through the early part of the transition to the nebular phase for completeness. 

a Angular size of the optical photosphere (i.e., the surface of last scattering). 
b Days since discovery, 1999 October 29 UT (HJD 2,451,480.94). 
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Table 8. EPM Distances to SN 1999em 



Method a 


Derived Distance 


Time of Explosion 13 




(Mpc) 


(d) 




7.68 ± 0.24 


-3.70 ±0.72 


BVI wcak 


8.27 ±0.18 


-5.69 ±0.50 


^-^weak 


8.79 ±0.28 


-6.51 ±0.77 


B V^trong 


8.18 ±0.41 


-3.83 ±1.15 


BV ^strong 


9.20 ±0.27 


-6.08 ± 0.74 


U/ s trong 


9.72 ± 0.42 


-7.16 ± 1.01 



Note. — The EPM distance derived to SN 1999em 
using the flux dilution factors of Eastman, Schmidt, & 
Kirshner (1996) as updated by Hamuy et al. (2001). 
Uncertainties are the la uncertainties resulting from 
1000 simulated sets of data characterized by the un- 
certainties in 9 and v found in Tables 6 and 7. 

a The broadband filter combination (BV, BVI, or 
VI) used to derive the photospheric temperature, and 
the absorption lines used to derive photospheric veloc- 
ity, where "strong" indicates Fe II AA4924, 5018, 5169 
and "weak" indicates the weaker, unblended features 
discussed in the text. 

b Days before discovery, 1999-10-29 UT (HJD 
2,451,480.94). 
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Fig. 20. — The parameters used to determine the distance to SN 1999em (equation [3]), with the 
line of best fit determined by using the criterion of least absolute deviations. The slope of this line 
yields the distance and the y-intercept the time of explosion of SN 1999em. Here, photospheric 
temperature is estimated using BV photometry, and photospheric velocity is determined using 
"Wweak; the velocity derived from weak, unblended features in the optical spectrum. 
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Fig. 21. — As in Figure 20, except using BVI photometry to determine photospheric color tem- 
perature. 
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Fig. 22. — As in Figure 20, except using VI photometry to determine photospheric color temper- 
ature. 



-62- 




Fig. 23. — As in Figure 20, except using t> s trong, the velocity determined by using the Fe II 
AA4924, 5018, 5169 line features, to estimate photospheric velocity. 
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Fig. 24. — As in Figure 23, except using VI photometry to determine photospheric color temper- 
ature. 
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Fig. 25. — As in Figure 23, except using VI photometry to determine photospheric color temper- 
ature. 
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found from the weak lines. Rounding to the nearest tenths place (in Mpc and days), then, yields 
D = 8.2 ± 0.6 Mpc with an explosion date of 5.3 ± 1.4 days prior to discovery as our best distance 
estimate. 

4. Discussion 

4.1. Sources of Systematic Uncertainty 

Our preferred EPM distance to SN 1999em is in good agreement with the recent distance 
estimate to NGC 1637 of D = 7.8^ Mpc found by Sohn & Davidge (1998) using the galaxy's 
brightest supergiant stars as well as the EPM distance of 7.5 ±0.5 Mpc recently reported by Hamuy 
et al. (2001) to SN 1999em. The major difference in the application of the EPM between our study 
and that of Hamuy et al. (2001) lies in the technique used to estimate photospheric velocity: 
while we rely on the velocity indicated by the flux minima of weak line features, Hamuy et al. 
cross-correlate the observed spectra with the model spectra of E96. A detailed discussion of the 
differences between these two methods is beyond the scope of the present study; we do note that the 
good agreement between the two resulting distances suggests that systematic differences between 
the two techniques may be small. It is also encouraging that two completely independent sets of 
data yield such consistent results. 

From the scatter among the EPM distances derived using the 3 different filter combinations, 
it is clear that systematic effects dominate the uncertainty in the distance to this well-observed 
SN, and are at least of order 10% — 20%. One obvious potential source of systematic error is the 
uncertainty in the flux dilution factor, £. The 5% — 10% uncertainty in ( (see § 1) reported for 
the range of supernova atmospheres modeled by E96 translates directly to distance uncertainties 
of 5% — 10%. The maximum deviation of any atmosphere studied by E96 from the mean £ is 
24%. If SN 1999em is such an extreme event and the degree of deviation between the true dilution 
factor and the average value varies among the different filter combinations studied, then the scatter 
among the derived distances may be explained. 

A second source of systematic error is the uncertainty in the extinction assumed for SN 1999em, 
which affects each bandpass differently. In general, EPM-derived distances are fairly robust to 
uncertainty in extinction (E96; Schmidt et al. 1992), due to the approximate cancellation of terms 
in equation (5): since dust both reddens and diminishes the light received, both the measured flux, 
f u , and the inferred blackbody flux, B U (T C ), decrease with increasing extinction. These two effects 
tend to cancel, producing a derived angular size, 6, that does not vary much with changes in the 
reddening. An additional point to consider, however, is that ( is also a strong function of T c for 
most filter combinations (generally decreasing with increasing T c ). The competition among these 
terms makes it difficult to generalize about how changes in reddening will alter the derived distance. 
For SN 1999em, lowering E(B—V) to 0.05 mag produces larger distances for all filter combinations 
{Dbv = 8.5 Mpc, Dbvi = 8.9 Mpc, and Dyi = 9.0 Mpc; all distances are derived using -u wea k f° r 
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the photospheric velocity [§ 2.4]), while raising the reddening to E(B— V) = 0.15 mag lowers the 
derived distance for the BV and BVI filter combinations {Dbv = 7.0 Mpc and Dbvi = 7.8 Mpc) 
and keeps Dyi at about the same distance {Dyi = 8.8 Mpc). For E(B— V) = 0.05 mag, taking 
the unweighted averages of the 3 filter combinations results in a distance and time of explosion 
of D = 8.8 ± 0.2 Mpc and t = 7.0 ± 1.3 days before discovery, while taking E(B-V) = 0.15 
mag yields D = 7.9 ± 0.9 Mpc and t a = 3.7 ± 1.6 days before discovery. For the large extinction 
(E[B—V] = 0.38) suggested by the strength of the interstellar Na I D lines (§ 2.3), we note that 
Dbv = 4.7 Mpc, Dbvi = 6.2 Mpc, and Dyi = 10.3 Mpc. The large discrepancies among the 
three distances in the high-extinction case, along with very poor fits obtained in the t vs. 9/v 
graphs (i.e., the points show large, systematic trends, and predict an explosion date somewhat 
after discovery) and unrealistically high temperatures derived at early times (Tbv > 10 6 K), all 
argue strongly against such a high reddening. While the lower extinction distances are much more 
self-consistent, we do note that no value of the extinction is found to make Dbv = Dbvi = Dyi- 

In addition to the uncertainty due to the dilution factor and extinction, asphericity could 
certainly contribute as well. However, it is not easy to see how asphericity could affect the individual 
bandpasses differently, to produce the filter-dependent discrepancies that are seen. L01 found that 
the polarization of SN 1999em increased with time, perhaps suggesting an increase in asphericity 
deeper into the ejecta. Barring a special viewing angle, increasing asphericity should produce 
a changing EPM distance with time. However, as shown in Figure 26, the derived distance to 
SN 1999em stays reasonably stable for all 3 filter combinations. Unless offsetting effects exist, the 
lack of obvious and consistent temporal trends across the three filter sets argues against the idea that 
the degree of asphericity increases with increasing depth in the SN ejecta. We therefore conclude 
that the most likely sources of systematic error are the uncertainty in the flux dilution factors 
and the extinction, and that these uncertainties clearly dominate over the statistical uncertainty 
resulting from the application of the EPM in any single filter combination alone. It should be noted 
that the uncertainty given for our preferred distance, D = 8.2 ± 0.6 Mpc, incorporates only the 
evident systematic uncertainty resulting from the application of the EPM with the three different 
filter combinations, and includes neither the uncertainty in extinction nor the disagreements among 
different groups on the overall level of the dilution factor (e.g., Baron et al. 1995). 



4.2. Comparison between EPM and Cepheid Distances 

Cepheid variables remain the most widely used and well studied extragalactic primary distance 
indicators. In particular, more than two dozen galaxies now have Cepheid distances obtained with 
the Hubble Space Telescope (HST) by the Distance Scale Key Project (e.g., Ferrarese et al. 2000; 
Freedman et al. 2001) and other groups (e.g., Saha et al. 1999; Tanvir, Ferguson, & Shanks 1999). 
Observing strategies and reduction techniques are well defined (e.g., Tanvir 1999). The slope and 
form of the Cepheid period-luminosity (PL) relation is quite well established from observations of 
LMC variables (e.g., Udalski et al. 1999). However, there is still some controversy over the distance 
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Fig. 26. — The derived EPM distance to SN 1999em as a function of time for the BV, BVI, and 
VI filter sets, with weak, unblended features used to estimate photospheric velocity. 
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to the LMC itself, which sets the PL zero-point (e.g., Silbermann et al. 1999). Considering all 
sources of error, the Key Project has recently announced "final" results for the Cepheid distances, 
with quoted random uncertainties of about 5% for each distance; systematic uncertainty, incor- 
porating uncertainty introduced by the LMC zero-point offset, metallicity correction, photometric 
zero-point, and aperture correction is at the 10% level (Freedman et al. 2001). 

Including the most recent additions, Cepheid distances now exist to 1 galaxy that hosted a 
SN II-P (SN 1973R), 5 galaxies that hosted photometrically or spectroscopically peculiar SNe II 
(SN 1970G, SN 1979C, SN 1987A, SN 1989L, and SN 1993J), and 3 galaxies in the same group 
in which a SN II-P occurred (SN 1968L, SN 1969L, and SN 1988A). Given the peculiar nature of 
some of these SNe and the uncertainties involved in comparing distances derived to different galaxies 
within the same group, it is quite impressive that E96 found -Dcepheids/-fEPM = 0.98 =L 0.08, based 
on the distances to 6 galaxies (or galaxies in the same group) that had Cepheid distances at the 
time of their study. Table 9 lists the Cepheid distances that were used by E96, the most recent 
Cepheid distances (taken from Freedman et al. 2001), and the EPM distances to 8 galaxies; we 
do not include SN 1993J in M81 since detailed models of its (likely unusual) early-time dilution 
factor evolution are needed to produce an EPM distance truely independent of the Cepheid scale 
(see Clocchiatti et al. 1995). 

Due to the recent modifications of the zero-point offset and refinement of the Cepheid distance 
technique itself, it is seen that the Cepheid distances to these galaxies have changed somewhat; 
in fact, nearly all of the Cepheid distances have decreased from their previous values. Using the 
revised Cepheid distances for the 6 SNe that were compared by E96, and using the same comparison 
technique (i.e., finding the simple weighted mean and la uncertainty in the weighted mean making 
no special consideration for the added uncertainty associated with comparing Cepheid distances 
to galaxies in the same group as a galaxy that hosted a SN II), now yields -Dcepheids/Z^EPM = 
0.88 ± 0.07; including all 8 available galaxies yields -Dcepheids/-CEPM = 0.87 ± 0.06. Of course, 
comparing group-member galaxy distances is a risky business and certainly increases the actual 
uncertainty of the comparison. If we eliminate the three group member comparisons (i.e., SN 1968L, 
SN 1969L, and SN 1988A), we then derive £> C c P hcids/£>EPM = 0.96 ± 0.09. Clearly, additional 
direct comparisons, particularly to spectroscopically normal SNe II-P, are needed to produce a 
more statistically meaningful comparison. Despite these being the "final" Cepheid distances, it 
should also be noted that other workers find evidence for a significantly longer Cepheid distance 
scale (e.g., Parodi et al. 2000), and that even within the Key Project itself longer distance scales 
have been given (e.g., Ferrarese et al. 2000). 

Nonetheless, the "best" Cepheid distances currently available suggest that EPM distances may 
be systematically too high. Since systematic uncertainty clearly remains at at least the 10% level for 
Cepheids, and at the 10 — 20% level for EPM, it is not clear how significant the observed discrepancy 
actually is. From our study of the effects that using strong lines to estimate photospheric velocity 
has on the derived distance, we might conclude that some of the difference between the EPM and 
Cepheid scales is due to inaccuracies in the photospheric velocity estimates; indeed, the distances 
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Table 9. EPM and Cepheid Distances to Type II Supernovae 



Supernova 


Galaxy Name 


Cepheid Distance 


EPM Distance 






SN host 


Cepheid galaxy 


Doid (Mpc) a 


D new (Mpc) b 


(Mpc) 


SN 


1968L 


NGC 5236 


NGC 5253 group 


4.1 ±0.4 


3.15 ±0.20 


4 -°-0.8 


SN 


1969L 


NGC 1058 


NGC 925 group 


9.29 ±0.69 


9.16 ±0.17 


10.6±i;? 


SN 


1970G 


NGC 5457 (M101) 


NGC 5457 (M101) 


7.5 ±0.8 


6.70 ± 0.34 




SN 


1973R 


NGC 3627 


NGC 3627 




10.05 ±0.37 


15 ±7 


SN 


1979C 


NGC 4321 (M100) 


NGC 4321 (M100) 


17.1 ± 1.8 


15.21 ±0.49 


15 ±4 


SN 


1987A 


LMC 


LMC 


0.051 ± 0.003 


0.050 ± 0.002 


0.049 ± 0.006 


SN 


1988A 


NGC 4579 


Virgo cluster 


17.1 ± 1.8 


15.28 ±0.35 


20 ±3 


SN 


1989L 


NGC 7331 


NGC 7331 




14.72 ±0.61 


17 ±4 



a Cepheid distance used by Eastman, Schmidt, & Kirshner (1996). 

b Metallicity-corrected Cepheid distance and random uncertainty reported by Freedman et al. (2001); 
systematic uncertainty in the Cepheid distance scale is ~ 10%. 

c Distances and uncertainties reported by Eastman, Schmidt, & Kirshner (1996). 
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derived by using the strong Fe II AA4924, 5018, 5169 lines would tend to produce EPM distances 
systmatically 5 — 10% too large. However, since it is not certain exactly which lines were used in 
each of the previous EPM studies to estimate photospheric velocity (and we note that using the 
Sc II A5658 line with A = 5657.9 A will produce velocities that are lower than those found from 
the weak, unblended lines [§ 2.2.3]), it is difficult to generalize. 

4.3. SNe II-P as Standard Candles 

To obtain the EPM distance to a SN II-P requires a well-spaced series of high-quality spectral 
and photometric data, and even with an excellent data set distance uncertainties of 10% — 20% 
currently still remain. As we look toward pushing EPM out to higher redshifts, it is clear that 
the most demanding observational requirement is the need to obtain multiple spectral epochs in 
order to determine the photospheric velocity. This is markedly different from the case for SNe la, 
which require only a single spectral epoch near (or, preferably, before) maximum light to identify 
it (e.g., Coil et al. 2000). An interesting suggestion has been made by Hofiich et al. (2000) that 
distances accurate to ~ 30% should be achievable by treating SNe II-P as standard candles. Based 
on theoretical models with a wide range of parameters (explosion energy, metallicity, mass loss 
of progenitor, etc.), Hofiich et al. (2000) find that although the peak luminosity varies greatly 
among SNe II-P, the mean absolute brightness during the plateau phase is quite insensitive to the 
initial parameters, with My (plateau) ps —17.6 ± 0.6 mag. Treating SNe II-P as standard candles 
is particularly attractive since they have a unique light-curve shape, obviating the need for even a 
single spectrum. Only a few deep images, roughly every 50 days (in the SN rest frame), are needed 
to discover, identify, and measure the mean plateau brightness; there is no need to follow the SN 
after the plateau ends and the SN becomes faint. 

To date, EPM distances have been derived to 10 SNe II-P, of which 8 have published F-band 
photometry sampling the plateau epoch. Adopting the EPM distance, the published photometry, 
and the extinction used in the EPM analysis (in some cases modified slightly due to the more 
accurate reddening values provided by the recent dust maps of Schlegel et al. [1998]), we derive the 
average F-band magnitude during the plateau for these 8 SNe II-P, and list the results in Table 10 
along with SN 1999em. To estimate the average F-band magnitude during the plateau, we took the 
average of the y-band data between days 20 and 100 after the date of explosion (as determined by 
the EPM analysis). With the exception of SN 1973R, which is anomalously bright, 23 all of the values 
cluster within 1.2 mag of each other. In fact, the weighted mean absolute F-band magnitude during 
the plateau epoch for SNe II-P, not including SN 1999em, is My (plateau) = — 16.8l ; 5 mag, where 
the reported uncertainty is the weighted standard deviation of the data from the mean. Including 



23 Note that the recent Cepheid distance to this galaxy is significantly less than that estimated by EPM (Table 9). 
If the Cepheid distance is correct, then SN 1973R actually has My (plateau) = —18.2 mag, making it somewhat more 
consistent with the other SNe II-P. 
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Table 10. Average y-Band Magnitude of SNe II-P During the Plateau Phase 



Supernova 






Distance 


M v c 


References 






I in a P" 1 


i TTIrlP' 1 


(MdcI 




Photometry Distance 


SN 1968L 


0.03 


12.0 


4.518:1 


-16.3±g;l 


1 


2 


SN 1969L 


0.18 


13.2 


10.611;? 


-16-9185 


1 


2 


SN 1973R 


2.80 


11.8 


15+? 


-19.1+J-8 

— l.U 


1 


2 


SN 1988A 


0.15 


14.6 


201s 


-16.918J 


1 


2 




1.0 


1 A 9 


18_ 2 


17 1+ - 2 


3 


2 


SN 1990ae 


0.5 


19.3 


11511 


-16.018;? 


1 


2 


SN 1992H 


0.28 


14.6 


20.2+j 


"16.9181 


4 


4 


SN 1992am 


0.3 


19.5 


180l| 


-16.818:1 


5 


2 


SN 1999em 


0.31 


13.7 


8 2+ - 6 


-15.918;! 


6 


6 



Note. — All uncertainties are la. 

a Estimated total extinction in the V band. If the reference from which the distance 
estimate was taken also included a derivation of the total extinction of the SN based 
on properties of the SN, then that value was used. In cases where the extinction was 
estimated for the Galaxy and host individually, and then summed, we have updated 
the Galactic extinction to be the values provided by the dust maps of Schlegel et al. 
(1998), replacing those given by Burstein & Heiles (1982). 

b Average apparent l/-band magnitude during the plateau phase, derived as the 
average of all F-band data between days 20 and 100 after explosion, corrected for 
extinction. 

c Average absolute y-band magnitude during the plateau phase, with uncertainty 
derived solely from the uncertainty in the EPM distance estimate. 

d Photometry data taken from Schmidt et al. (1992), and references therein. 

References. — (1) Schmidt et al. 1992, and references therein; (2) Schmidt et al. 
1994a; (3) Schmidt et al. 1993; (4) Clocchiatti et al. 1996; (5) Schmidt et al. 1994b; 
(6) This work. 
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SN 1999em lowers the average to My (plateau) = — 16.4lg g mag. The data thus favor a somewhat 
fainter mean plateau V-band brightness than the models of Hoflich et al. (2000) predict, but with 
a similar spread around the mean, much of which may be due to the uncertainty in the EPM 
technique itself. From just these 9 SNe II-P, we would conclude that distances to SNe II-P accurate 
to ±28% (la) could be achieved with the standard-candle assumption, which compares well with 
the current level of systematic uncertainty in EPM, and is much less observationally taxing. Clearly, 
additional distances to SNe II-P are needed to better quantify the uncertainty resulting from the 
standard-candle assumption. Although the use of SNe II-P as standard candles will not achieve the 
same individual accuracy as SNe la, the removal of the need to obtain a spectrum, coupled with 
their certain existence at high redshift (before SNe la even occur, perhaps), may yet make them 
the best SN class for cosmology at redshifts beyond 2.0. 



4.4. Is SN 1999em a "Typical" SN II-P? 

In § 2.3 we raised the possibility that the B — V color evolution of SN 1999em might be 
somewhat different from the normal evolution of SNe II-P during the plateau phase. In § 4.3 we 
found that SN 1999em was somewhat fainter during the plateau than other SNe II-P with EPM 
distances; at our derived distance, SN 1999em also reached a peak extinction-corrected i?-band 
magnitude of only Mb = —16.2 ± 0.2 mag, which places it on the fainter end of the range observed 
for SNe II-P (Patat et al. 1994; Miller & Branch 1990). This naturally leads to the question: is 
SN 1999em unusual? If so, could this have affected the EPM analysis? More to the point, are 
the observed properties of SN 1999em similar to the theoretical models used to derive the dilution 
factors (C) by E96? 

Since SN 1999em evolved more rapidly toward red colors than the prototypical Type II-P 
SN 1969L (see Figure 13), we first consider the temporal evolution of its color temperature. The 
lowest color temperatures for any of the models studied by E96 are those for model sl5.60.1, which 
has T B v = 4132 K and T B vi = 5418 K, and 10h.60.1, which has T VI = 6638 K; 24 only 1 other 
model studied (p6.40.1) has Tbv < 5000 K and Tbvi < 6000 K. However, a quick look at Figure 5 
and Table 7 reveals that SN 1999em drops below T BV = 5000 K and T BV i = 6000 K near day 
30, and below Tbv = 4132 K and Tbvi = 5418 K around day 50. This means that after about 
day 30, the color temperature of SN 1999em is only represented by 2 of the 63 models studied by 
E96, and after day 50 it has no representation in the model set used to derive C(T C ), but rather 
relies on extrapolation. At these low temperatures, (bv an d (bvi are changing very rapidly with 
temperature. A contribution to the distance uncertainty could thus come from the lack of models 
near the observed color temperature of SN 1999em during the late portion of the recombination 



24 We note that a small change in the model's color-temperatures will result from the use of the Hamuy et al. (2001) 
filter functions in place of those of E96; for the low-temperature regime of interest here, however, these differences 
should be minimal. 
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phase. We note that increasing the reddening of SN 1999em to the upper limit given by Baron et 
al. (2000; see also § 2.3) of E(B-V) = 0.15 mag raises T BVI by 400 K and T BV by only about 
200 K on day 30. 

From its relative faintness and rapid color evolution during the plateau there is therefore some 
evidence that SN 1999em was a somewhat unusual II-P event. It is possible that the average £ 
values derived by E96 do not adequately represent the evolution seen in SN 1999em, and may point 
toward the need for models specifically crafted for it to increase the accuracy of its distance. 

One way to help place the theoretically derived flux dilution factor on a more sound empirical 
footing is to directly compare the EPM distance to NGC 1637 with that derived using other primary 
extragalactic distance indicators. At a distance of 8.2 Mpc, NGC 1637 makes an ideal target for 
a Cepheid study using HST: it is rather face-on (i « 32°), has low foreground extinction, and has 
a metallicity of log (O/H) + 12 ~ 9.1 in the region near the SN (van Zee et al. 1998), typical of 
other HST galaxies for which Cepheid distances have been derived (e.g., Ferrarese et al. 2000; 
see also Tanvir 1999). An approved Cycle 10 HST program to derive a Cepheid-based distance to 
NGC 1637 is currently underway, and the results will be directly compared with the EPM distances 
derived to SN 1999em. 



5. Conclusions 

We present 30 optical spectra and 49 photometric epochs of SN 1999em sampling the first 517 
days of its development. SN 1999em displays the major photometric and spectral characteristics 
of a SN II-P, exhibiting a distinct photometric plateau lasting about 100 days from explosion, and 
well-developed P-Cygni spectral line features. We suspect that SN 1999em suffers from minimal 
reddening, and concur with the range of possible values found by Baron et al. (2000), and adopt 
E(B-V) = 0.10 ± 0.05 mag. 

We find a number of interesting features in the spectral evolution of the P-Cygni line profiles, 
particularly Ha. We identify the dominant ions responsible for most of the absorption features seen 
in the optical portion of the spectrum of SN 1999em during the plateau phase. The photospheric 
velocity derived from four weak, unblended absorption features is found to be 5% — 10% lower 
than that derived using the stronger Fe II AA4924, 5018, 5169 lines, which have traditionally 
been used to estimate the photospheric velocity in EPM studies. The velocity offset between the 
weakest lines and the stronger Fe II features is shown to exist in other SNe II as well. If this is a 
generic feature of SNe II and weaker lines are indeed a better tracer of photospheric velocity, then 
EPM estimates previously made with the stronger Fe II AA4924, 5018, 5169 lines are skewed by 
5% — 10% toward larger distances. We note that the comparison of EPM distances with the most 
recent Cepheid distances yields -Dcepheids/-DEPM = 0.87 ± 0.06, a discrepancy that is in the same 
sense (i.e., -Depm > -Dcepheids) as the aforementioned velocity offset; difficulties in this comparison 
exist, however, due to the peculiar nature of several of the SNe II, the use of Cepheid distances to 
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galaxies that are only in the same group in which the SN occurred (i.e., not the host galaxy itself), 
and existing discrepancies in the Cepheid distance scale. 

We derive an EPM distance to SN 1999em using the flux dilution factors (() of E96 as modified 
by Hamuy et al. (2001). We use three different filter combinations (BV, BVI, and VI) to determine 
the photospheric color temperature. The total spread in the distances derived using the different 
filter combinations is about 15% (when using the weak, unblended absorption features to estimate 
photospheric velocity), in the sense Dyi > Dbvi > Dbv- These discrepancies cannot be explained 
by statistical uncertainty in the measured quantities, and must be systematic in nature, most likely 
due to uncertainty in the dilution factor and perhaps aggravated by the lack of models used to derive 
the dilution factor with properties similar to SN 1999em: SN 1999em became redder more quickly 
than the prototypical Type II-P SN 1969L, and reached color temperatures lower than any model 
studied by E96 about half way through the plateau phase. Taking the observed scatter among 
the different filter combinations into account, our best estimate for the distance to SN 1999em is 
D = 8.2±° 6 6 Mpc, with an explosion date of t Q = HJD 2,451,475.6^;^, or 5.3 ± 1.4 days before 
discovery. This leads to an extinction-corrected peak 5-band magnitude of Mb = —16.2 ±0.2 mag, 
and an average U-band plateau brightness of My(plateau) = —15.9 ± 0.2 mag, which is somewhat 
fainter than average SNe II-P. 

Finally, we investigate the theoretical prediction by Hdflich et al. (2000) that SNe II-P should 
exhibit a rather small dispersion in their average U-band magnitude during the plateau phase. 
Their unique light-curve shape (making spectral identification unnecessary) and certain existence 
at high redshift make them attractive cosmological tools. From 8 SNe II-P with previously published 
distances and SN 1999em, we find My (plateau) = — 16.4^'g. If this small sample is representative 
of all SNe II-P, then distances to SNe II-P accurate to ~ 30% (Id) may be possible without having 
to acquire the extensive data set required for an EPM analysis, thereby increasing the viability of 
this class of objects as cosmological beacons at z > 2. 
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